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I. INTRODUCTION 

Cosmological inflation [1] has become the dominant paradigm to understand the initial conditions for the CMB 
anisotropies and structure formation. In the inflationary picture, the primordial cosmological perturbations are 
created from quantum fluctuations "redshifted" out of the horizon during an early period of accelerated expansion of 
the universe, where they remain "frozen". They are observable as temperature anisotropies in the CMB. This picture 
has recently received further spectacular confirmation by the Wilkinson Microwave Anisotropy Probe (WMAP) three 
year set of data [2]. Since the observed cosmological perturbations are of the order of 10 -5 , one might think that 
first-order perturbation theory will be adequate for all comparisons with observations. This might not be the case, 
though. Present [2] and future [3] experiments may be sensitive to the non-linearities of the cosmological perturbations 
at the level of second- or higher-order perturbation theory. The detection of these non-linearities through the non- 
Gaussianity (NG) in the CMB [4] has become one of the primary experimental targets. 

A possible source of NG could be primordial in origin, being specific to a particular mechanism for the generation 
of the cosmological perturbations. This is what makes a positive detection of NG so relevant: it might help in 
discriminating among competing scenarios which otherwise might be undistinguishable. Indeed, various models of 
inflation, firmly rooted in modern particle physics theory, predict a significant amount of primordial NG generated 
either during or immediately after inflation when the comoving curvature perturbation becomes constant on super- 
horizon scales [4]. While single-field [5] and two(multi)-field [6] models of inflation generically predict a tiny level of 
NG, 'curvaton-type models', in which a significant contribution to the curvature perturbation is generated after the 
end of slow-roll inflation by the perturbation in a field which has a negligible effect on inflation, may predict a high 
level of NG [7, 8] . Alternatives to the curvaton model are those models characterized by the curvature perturbation 
being generated by an inhomogeneity in the decay rate [9, 10], the mass [11] or the interaction rate [12] of the particles 
responsible for the reheating after inflation. In that case the reheating can be the first one (caused by the scalar field(s) 
responsible for the energy density during inflation) or alternatively the particle species causing the reheating can be a 
fermion [13]. Other opportunities for generating the curvature perturbation occur at the end of inflation [14], during 
preheating [15], and at a phase transition producing cosmic strings [16]. 
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On the other hand there exist many sources of NG in the CMB anisotropics beyond the primordial ones, which is 
essential to characterize in order to distinguish them from a possible primordial signal. For example, statistics like the 
bispectrum and the trispectrum of the CMB can be used to assess the level of primordial NG on various cosmological 
scales and to discriminate it from the one induced by secondary anisotropies and systematic effects [4, 17, 18, 19]. In 
this case a positive detection of a primordial NG in the CMB at some level might therefore confirm and/or rule out 
a whole class of mechanisms by which the cosmological perturbations have been generated. 

Therefore it is of fundamental importance to provide accurate theoretical predictions of all the possible non-linear 
effects contributing to the overall NG in the CMB anisotropies. At second-order in the perturbation theory one 
should provide a full prediction for the second-order radiation transfer function. A first step towards this goal has 
been taken in Ref. [20] (see also [21]) where the full second-order radiation transfer function for the CMB anisotropies 
on large angular scales in a flat universe filled with matter and cosmological constant was computed, including the 
second-order generalization of the Sachs- Wolfe effect, both the early and late Integrated Sachs- Wolfe (ISW) effects 
and the contribution of the second-order tensor modes. These effects are due to gravity. In Ref. [22] we presented the 
computation of the full system of Boltzmann equations at second-order describing the evolution of the photon, baryon 
and CDM fluids, neglecting polarization. In this way we accounted also for the small scale effects due to the collision 
terms. The equations we derived allow to follow the time evolution of the CMB anisotropies at second-order at all 
angular scales from the early epoch, when the cosmological perturbations were generated, to the present through the 
recombination era. Ref. [22] sets therefore the stage for the computation of the full second-order radiation transfer 
function at all scales and for a a generic set of initial conditions specifying the level of primordial non-Gaussianity. 

At second-order one can see that there are many sources of NG in the CMB anisotropies, beyond the primordial one. 
The most relevant sources are the so-called secondary anisotropies, which arise after the last scattering epoch. The 
so called scattering secondaries include the thermal Sunyaev-Zel'dovich effect, where hot electrons in clusters transfer 
energy to the CMB photons, the kinetic Sunyaev-Zel'dovich effect produced by the bulk motion of the electrons 
in clusters, the Ostriker-Vishniac effect, produced by bulk motions modulated by linear density perturbations, and 
effects due to reionization processes. Gravitational secondaries are effects mediated by gravity and include the change 
in energy of photons when the gravitational potential is time-dependent or the gravitational lensing. Secondaries 
that result from a time-dependent potential arc the ISW produced mainly on large scales when the dark energy 
at late times becomes dominant and the potential starts to decay, or the Rees-Sciama effect, produced during the 
matter-dominated epoch at second-order and by the time evolution of the potential on non- linear scales. Gravitational 
lensing which causes the deflection of the photons'path from the last scattering to us, does not create anisotropies, 
but it only modifies existing ones. All of these secondaries effects are most significant on small angular scales (except 
for the ISW effect). Moreover the three-point function arising from the correlation of the gravitational lensing effect 
and the ISW effect generated by the matter distribution along the line of sight [23, 24] and the Sunyaev-Zel'dovich 
effect [25] arc large and detectable by Planck [26]. Of course on small angular scales, fully non- linear calculations 
of specific effects like Sunyaev-Zel'dovich, gravitational lensing, etc. would provide a more accurate estimate of the 
resulting CMB anisotropy, however, as long as the leading contribution to second-order statistics like the bispectrum 
is concerned, second-order perturbation theory suffices. 

In this paper we will focus on another relevant source of NG: the non-linear effects operating at the recombination 
epoch. The dynamics at recombination is quite involved because all the non-linearities in the evolution of the baryon- 
photon fluid at recombination and the ones coming from general relativity should be accounted for. The present paper 
can be considered as an application of the equations found in Ref. [22] and it offers an analytical study at second-order 
of this complicated dynamics. This allows to account for those effects that at the last scattering surface produce a 
non-Gaussian contribution to the CMB anisotropies that add to the primordial one. Such a contribution is so relevant 
because it represents a major part of the second-order radiation transfer function which must be determined in order 
to have a complete control of both the primordial and non-primordial part of NG in the CMB anisotropies and to 
gain from the theoretical side the same level of precision that could be reached experimentally in the near future [4]. 

In order to achieve this goal, we have considered the Boltzmann equations derived in Ref. [22] at second-order 
describing the evolution of the photon, baryon and CDM fluids, and we have manipulated them further under the 
assumption of tight coupling beteween photons and baryons. This leads to the generalization at second-order of 
the equations for the photon energy density and velocity perturbations which govern the acoustic oscillations of the 
photon-baryon fluid for modes that are inside the horizon at recombination. The evolution is that of a damped 
harmonic oscillator, with a source term which is given by the gravitational potentials generated by the different 
species. An interesting result is that, unlike the linear case, at second-order the quadrupole moment of the photons 
is not suppressed in the tight coupling limit and it must be taken into account. We also find that the second-order 
CMB anisotropies generated at last scattering do not reduce only to the energy density and velocity perturbations of 
the photons evaluated at recombination, but a number of second-order corrections due to gravity at last scattering 
arise from the Boltzmann equations of Ref. [22]. We compute them when we decompose the CMB anisotropies in 
multipole moments. 
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The bulk of the paper deals with the computation of the analytical solutions for the acoustic oscillations of the 
photon-baryon fluid at second-order. These solutions arc derived adopting some simplifications which are also standard 
for an analytical treatment of the linear CMB anisotropies, and which nonetheless allow to catch most of the physics at 
recombination. One of these simplifications is to study separately two limiting regimes: intermediate scales which enter 
the horizon in between the equality epoch {rj eq ) and the recombination epoch (r] r ), with i]^ 1 <C k <C r\^ q x , and shortwave 
perturbations, with k 3> rj~ q ~, which enter the horizon before the equality epoch. An alternative approach could be to 
derive a semianalytical solution by using the fits of Ref. [28] for the linear gravitational potentials. Otherwise a full 
numerical evaluation can be performed [29] using the set of Boltzmann equations of Ref. [22]. However in this paper 
our main concern is to provide a simple estimate of the quantitative behaviour of the non-linear evolution taking 
place at recombination, offering at the same time all the tools for a more accurate computation. Notice that the case 
k 3> rj^ 1 has been treated in two steps. First we just assume a radiation dominated universe, and then we give a 
better analytical solution by solving the evolution of the perturbations from the equality epoch onwards taking into 
account that the dark matter perturbations around the equality epoch tend to dominate the second-order gravitational 
potentials. As a byproduct, this last step provides the Meszaros effect at second-order. In deriving the ananlytical 
solutions we have accurately accounted for the initial conditions set on superhorizon scales by the primordial non- 
Gaussianity. In fact the primordial contribution is always transferred linearly, while the real new contribution to the 
radiation transfer function is given by all the additional terms provided in the source functions of the equations. Let 
us stress here that the analysis of the CMB bispectrum performed so far, as for example in Ref. [2, 26, 30], adopt 
just the linear radiation transfer function (unless the bispectrum originated by specific secondary effects, such as 
Rees-Sciama or Sunyaev-Zel'dovich effects, are considered). 

Since the paper serves for different purposes and achieves different goals, we summarize them in the following: 

• We compute the second-order CMB anisotropies generated by non-linearities at recombination which will add 
to the primordial non-Gaussianity. 

• We provide analytical solutions for acoustic oscillations of the photon-baryon fluid at second-order in the tight 
coupling limit starting from the Boltzmann equations derived in Ref. [22]. 

• We compute the evolution of the CDM density perturbations (and the gravitational potentials) accounting for 
those modes that enter the horizon during the radiation dominated epoch. This allows to determine the second- 
order transfer function for the density perturbations, and in particular the generalization at second-order of the 
Meszaros effect. 

• We provide the multipolc moments of the CMB anisotropies and hence we are able to compute that part 
of the second-order radiation transfer function that corresponds to small-scale effects at recombination, thus 
complementing the results of Ref. [20] . 

The paper is organized as follows. In Section II we just report the Boltzmann equations for the photons derived in 
Ref. [22]. In Sec. Ill wc recall how to treat them in the tight coupling limit at linear order in the perturbations and 
the standard way to get the analytical solutions for the photon-baryon fluid. In Section IV we derive the equations 
for the second-order energy density and velocity perturbations of photons in the tight coupling limit. A subsection is 
devoted to compute the second-order quadrupole moment of the photons. Sec. V deals with the expansion of CMB 
anisotropies in multipole moments and with the computation of those contributions to the CMB anisotropies which 
are generated at recombination. In Sec. VI we derive the analytical solutions describing the acoustic oscillations at 
second-order of the photon baryon fluid. Wc derive these solutions accounting for the primordial non-Gaussianity, and 
in the two regimes r/^ 1 <C k <C r/^ 1 (sec. VII) and k 3> r/^ 1 (Sec. VIII). In Sec. IX we compute at first- and second- 
order the evolution on subhorizon scales of the density perturbations of CDM, thus arriving at a generalization of the 
Meszaros effect. This result also allows to give a refined prediction for the CMB anisotropies in the case k ^> f]~ q . Sec. 
X contains our conclusions, and we also provide some Appendices which mianly treat the gravity perturbations, and 
where we provide the generic solutions for the evolution of the second-order gravitational potentials for a radiation- 
or matter dominated universe. 



II. THE BOLTZMANN EQUATIONS 

In this Section we just report the Boltzmann equations derived in Ref. [22], while the goal of Sec. Ill and IV is 
to derive the moments of the Boltzmann equations for photons in the limit when the photons are tightly coupled to 
the the baryons (the electron-proton system) due to Compton scattering. We will first review briefly the standard 
computation at linear order and then derive the equations at second-order in the perturbations, pointing out some 
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interesting differences with respect to the linear case. The starting point is the Boltzmann equation at first- and 
second-order [22] 
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Let us recall some definitions of the quantities appearing in Eqs. (2.1)-(2.2). $ = $W + $( 2 )/2 and * = *W + 
\I>( 2 )/2 are the gravitational potentials in the Poisson gauge, while Wj and Xij are the second-order vector and tensor 
perturbations of the metric according to Eq. (A.l). The photon temperature anisotropies are given by 
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which represents the photon fractional energy perturbation (in a given direction) being the integral of the photon 
distribution function / = + /( 2 ) /2 over the photon momentum magnitude p (p l = pn l ). The angular dependence 
of the photon anisotropies A can be expanded as 
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where we warn the reader that the superscript stands by the order of the perturbation, while the subscripts indicate 
the order of the multipoles. At first order one can drop the dependence on m setting m = so that a£^ = 

(—i)~ e (2£ + l)S m0 A^\ It is understood that on the left-hand side of Eq. (2.2) one has to pick up for the total time 
derivatives only those terms which contribute to second-order. Thus we have to take (see Ref. [22]) 
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Notice that we can write dA^/dn 1 = {dA^ / dx l ){dx l / dn 1 ) = (dA^ /dx l )(n - n t ). 

In Eq. (2.2) 6^ is the relative energy density perturbation of the electrons. These are in turn strongly coupled 
with protons (p) via Coulomb interactions, such that the density constrasts and the velocities are driven to a common 
value S e = <L 



up = 8b and v e = v& = v for what can then be called the baryon fluid. Finally 



(2.7) 



is the differential optical depth for the Compton scatterings between photons and free electrons, with gt the Thomson 
cross section, a the scale factor, and n e represents the mean density of free electrons. The tightly coupled limit 
corresponds to the Compton interaction rate much bigger than the expansion of the universe, t' /H » 1 (or t > 1). 
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III. LINEAR BOLTZMANN EQUATIONS 

The first two moments of the photon Boltzmann equations are obtained by integrating Eq. (2.1) over df2 n /47r and 
dQ n n l I A-k respectively and they lead to the density and velocity continuity equations 

Aoo' + ^f -#«' = 0, (3.1) 



v^' + ^n«^ + La®-' + *m = r> (V 1 )* v w) . (3.2) 

Here we recall that 5^ = Aqq = J dilA^/An and that the photon velocity is given by [22] 
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IP- 7 is the quadrupole moment of the photons defined as 
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The two equations above are complemented by the momentum continuity equation for baryons, which can be conve- 
niently written as 



T 

where we have introduced the baryon-photon ratio 



(3.5) 



i?=|^. (3.6) 

Equation (3.5) is in a form ready for a consistent expansion in the small quantity r _1 which can be performed in the 

tight coupling limit. By first taking v^ z = v^ 1 at zero order and then using this relation in the L.H.S. of Eq. (3.5) 
one obtains 



(3.7) 



Such an expression for the difference of velocities can be used in Eq. (3.2) to give the evolution equation for the photon 
velocity in the tightly coupled limit 

v^' + n-^-v^ + + d>w = o . (3.8) 

7 l + i? 7 41 + i? 

Notice that in Eq. (3.8) we are neglecting the quadrupole of the photon distribution n' 1 )*- 7 (and all the higher moments) 
since it is well known that at linear order such moment (s) are suppressed in the tight coupling limit by (successive 
powers of) 1/r with respect to the first two moments, the photon energy density and velocity. Eqs. (3.1) and (3.8) are 
the master equations which govern the photon-baryon fluid acoustic oscillations before the epoch of recombination 
when photons and baryons are tightly coupled by Compton scattering. 

In fact one can combine these two equations to get a single second-order differential equation for the photon energy 
density perturbations Aqq . Deriving Eq. (3.1) with respect to conformal time and using Eq. (3.8) to replace diV^ 1 
yields 

(A«" - 4*«") + H-^ (Ag* - 4*W') - c 2 V 2 (A« - 4*«) = lv 2 (V> + ^) , (3.9) 

where we have introduced the photon-baryon fluid sound of speed c s = 1/ ^/3(1 + R). In fact in order to solve 
Eq. (3.9) one needs to know the evolution of the gravitational potentials. We will come back later to the discussion 
of the solution of Eq. (3.9). 
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A useful relation we will use in the following is obtained by considering the continuity equation for the baryon 
density perturbation. By perturbing at first-order Eq. (6.22) of Ref. [22] we obtain 

tfp + vji - 3* W = . (3.10) 

Subtracting Eq. (3.10) form Eq. (3.1) brings 

Aff'-J^ + l^-^-O, (3.11) 
which implies that at lowest order in the tight coupling approximation 

Aft = f , (3.12) 

for adiabatic perturbations. 



A. Tightly coupled solutions for linear perturbations 

In this section we briefly recall how to obtain at linear order the solutions of the Boltzmann equations (3.9). These 
correspond to the acoustic oscillations of the photon-baryon fluid for modes which are within the horizon at the time 
of recombination. It is well known that, in the variable (Aqq — 4\E'( 1 )), the solution can be written as [28, 31] 

[1 + i?(r 7 )] 1 / 4 (AW - 4*W) = Acos[fcr s (r/)] + B sin[fcr,(»y)] 

- 4 A |" ^'[1 + fl(„')] 3 / 4 ^'(r,') + sinlK^,,) - r s („'))] (3.13) 

where the sound horizon is given by 

r s ( V ) = ["dr/c a (r/), (3.14) 
Jo 

with the ratio R defined in Eq. (3.6). The first line of Eq. (3.13) corresponds to the solutions of the homogeneous 
equation, while the remaining integral corresponds to a particular solution of Eq. (3.13). The constants A and B 
must be fixed according to the initial conditions. 

In order to give an analytical solution that catches most of the physics underlying Eq. (3.13) and which remains 
at the same time very simple to treat, we will make some simplifications following Ref. [32, 33]. First, for simplicity, 
we are going to neglect the ratio R wherever it appears, except in the arguments of the varying cosines and sines, 
where we will treat R — R* as a constant evaluated at the time of recombination. In this way we keep track of a 
damping of the photon velocity amplitude with respect to the case R = which prevents the acoustic peaks in the 
power-spectrum to disappear. Treating R as a constant is justified by the fact that for modes within the horizon the 
time scale of the ocillations is much shorter than the time scale on which R varies. If R is a constant the sound speed 
is just a constant 

(3.15) 



S v/3(l + i?,) ' 

and the sound horizon is simply r s (rf) — c s rj. 

Second, we are going to solve for the evolutions of the perturbations in two well distinguished limiting regimes. One 
regime is for those perturbations which enter the Hubble radius when matter is the dominant component, that is at 
times much bigger than the equality epoch with k -C k eq ~ rf^ q x , where k eq is the wavenumber of the Hubble radius at 
the equality epoch. The other regime is for those perturbations with much smaller wavelenghts which enter the Hubble 
radius when the universe is still radiation dominated, that is perturbations with wavenumbers k 3> k eq ~ r/~ q . In fact 
we are interested in perturbation modes which are within the horizon by the time of recombination r\* . Therefore we 
will further suppose that 77* » rj eq in order to study such modes in the first regime. Even tough 77, » ij eq is not the 
real case, it allows to give some analytically solutions. 

Before solving for these two regimes let us fix the initial conditions which are taken on large scales deep in the 
radiation dominated era (for 77 — > 0). During this epoch, for adiabatic perturbations, the gravitational potentials 
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remain constant on large scales (we are negelcting anisotropic stresses so that ~ Vl/W) and from the (0 — 0)- 
componcnt of Einstein equations 

^(0) = - 1 -Ag ) (0). (3.16) 

On the other hand, from the energy continuity equation (3.1) on large scales 

A$ - 4* (1) = const. , (3.17) 

from which the constant is fixed to be const. = — 6$^^(0) and thus we find B = and A = -6$«(0). 
With our semplifications Eq. (3.13) then reads 

(A$ - 4* (1) ) = -6$ (1) (0)cos(cj ?7) - 8-^= C H^{rj[) sin[w (»7 - v')} , (3-18) 

v3 Jo 

where luq = kc s . 



B. Perturbation modes with k <C k eq 

This regime corresponds to perturbation modes which enter the Hubble radius when the universe is matter domi- 
nated at times r\ 3> rj eq . During matter domination the gravitational potential remains constant (both on super- horizon 
and sub-horizon scales), as one can see for example from Eq. (B.l), and its value is fixed to $^^(k,ry) = ^$^^(0), 
where $^^(0) corresponds to the gravitational potential on large scales during the radiation dominated epoch. Since 
we are interested in the photon anisotropics around the time of recombination, when matter is dominating, we can 
perform the integral appearing in Eq. (3.13) by taking the gravitational potential equal to its value during matter 
domination so that it is easily computed 

2 / eV$ (1 V)sinM*7-r/)] = ^ ^ (1 - cos(^)) . (3.19) 

Thus Eq. (3.18) gives 

A$ -4*W = ^$ (1) (0)cos(c^) - ^$«(0). (3.20) 

The baryon-photon fluid velocity can then be obtained as d l v^ 1)l = -3(A$ - 4*W)'/4 from Eq. (3.1). In Fourier 
space 

ih v^ 1 = (0) sin(w r/)wo , (3.21) 

where we use the convention diV^ % — > iki ^^(k) or equivalentely 

< }i = A$(i)( ) M«X>v)c. , (3.22) 

since the linear velocity is irrotational. 



C. Perturbation modes with k » k eq 

This regime corresponds to perturbation modes which enter the Hubble radius when the universe is still radiation 
dominated at times rj <c r) eq . In this case an approximate analytical solution for the evolution of the perturbations 
can be obtained by considering the gravitational potential for a pure radiation dominated epoch, given by Eq. (B.16). 
For the integral in Eq. (3.18) we thus find 

f $ (1) (t/) sinM?7 - rf)\ = -^-cos(w oV ) , (3.23) 
Jo Zu} 
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where we have kept only the dominant contribution oscillating in time, while neglecting terms which decay in time. 
The solution (3.18) becomes 

A$ - 4*W = 6$ (1) (0) cos(w r?) , (3.24) 

and the velocity is given by 

< )4 = -4^ (1) (0)sm(^)c s , (3.25) 

Notice that the solutions (3.24)-(3.25) arc in fact correct only when radiation is dominating. Indeed between 
the epoch of equality and recombination, matter will start to dominate. We will account for such a period and its 
consequences on the CMB anisotopy evolution in a separate section showing that some corrections must be properly 
taken into account. However for the time being we will keep on discussing the case k ^> k eq just by adopting the 
gravitational potential for a radiation dominated epoch, since it can be considered a first useful approximation in 
order to give the main quantitative features. 

Before moving to the study of the tightly coupled solutions for the second-order CMB anisotropies, we want 
to recover the solutions (3.24)-(3.25) in an alternative way, which will be particularly useful for the second-order 

case. Instead of solving the second-order differential equation (3.9) for Aqq\ we use directly the energy continuity 
equation (3.1). The reason is that for the case of radiation domination the gravitational potential (B.16) at late times 
decays as rf~ 2 being approximated by 

Notice that in the expression for the gravitational potential (B.16) we account for the sound-speed of the photon- 
baryon fluid, and as usual we keep it only in the argument of the sines and cosines. 

On the other hand from the (0 — z)-component of Einstein equation (B.17) we find that 

<>' * -^d^ (1V = -4* (1)(0) T s[n ^oV)cs , (3.27) 



and its divergence 



8 % v^ ~ - JLv 2 *' 1 '' = ^«(0)fcsin(^K , (3.28) 



where the second equalities are written in Fourier space and we have kept only the dominant terms at late time scaling 
like sin(wofy). Notice that we recover the same result of Eq. (3.25). As a result in Eq. (3.1) the gravitational potential 
can be neglected so that 

C-^^VV)'. (3.29) 
We integrate Eq. (3.29) using the late time expression (3.26) for the gravitational potential to find 

A$ = 6$ (1) (0)cos(w 77) . (3.30) 

The result in Eq. (3.30) agrees with the previous result (3.24) since the gravitational potential can be neglected at 
late times. 



IV. SECOND-ORDER BOLTZMANN EQUATIONS IN THE TIGHTLY COUPLED LIMIT 

Let us now treat in a similar way the photon Boltzmann equations at second-order in the cosmological perturbations 
exploiting the regime of tight coupling between the photons and the baryons to find the governing equations for the 
acoustic oscillations of the photon-baryon fluid at second-order. While we already know that the L.H.S. of the 
equations at second-order will have the same form as for the linear case, one of the main points here is to compute the 
source term on the R.H.S. of the moments of the Boltzmann equations which will consist of fist-order squared terms. 
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A. Energy continuity equation 

Let us start by integrating Eq. (2.2) over (Kln/in to get the evolution equation for the second-order photon energy 
density perturbations A$ 

Aqq" + IdivW + ^ (a«^) 4*( 2 )' + + ^)d iV ^ + 2(r? - % )($ (1) + * (1) ),, s i n( 1 )« 

+ ^VA&fo - r,) - 8*«'A« + f = - |/t,« - , (4.1) 

where we have used the explicit dchnition for the second-order velocity of the photons [22] 



4 v) 



3 2 



^A<»>„<-!W. (4.2) 



We can now make use of the tight coupling expansion to simplify Eq. (4.1). In the L.H.S. we use diV^ 1 — diV^ 1 — 
3^(1)' _ flW' = 3^(1)' _ 3A^ ); /4 obtained in the tightly coupled limit from Eqs. (3.10) and (3.12). On the other hand 
in the R.H.S. of Eq. (4.1) one can write 

(v™ - vW) = £ (vW + + *™) = f (-^<>* - , (4-3) 

by using Eq. (3.7) and the evolution equation for the photon veleocity (3.8). We thus arrive at the following equation 

A^' + ld lV ^-4^y=S A , (4.4) 

where the source term is given by 

S A = (A« 2 )' - + ^>)(4*«' - A«') - §t,W'(A&> + 4$«), 4 + |(„ - + 

_ |i2 ( -^-^) 2 - i^^^) . (4.5) 



B. Velocity continuity equation 

We now derive the second moment of the Boltzmann equation (2.2) and then we take its tight coupling limit. The 
integration of Eq. (2.2) over dtt n yields the continuity equation for the photon velocity 

A $(i),i($d) + vp(D) _ ($(D + *(i))a,.naw* - + *( 1 ))a« ,< + ^(r? - ^)(* (1) + * (1) ) ^'a,-^ 



-fa- 



+2<;«n( 1)j 



(4.6) 



The difference between the second-order baryon and photon velocities (t;( 2 ) 1 — v^ 1 ) appearing in Eq. (4.6) is obtained 
from the baryon continuity equation which can be written as (see Ref. [22]) 



„(2)i/ + Hv (2)i + 2(J i, + m ^ + $ (2),i^ _ 2$ (l)/ w (l)i + ^(1)2 + 2$ (l),i ($ (l) + ^(D)] _ ^(l) n (l)ii 

2(AW+$W)( 1) ( 1 )'- B f). (4.7) 



v (2)i = v (2), + ^ 

' r' 
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We want now to reduce Eq. (4.6) in the tightly coupled limit. We first insert the expression (4.7) in Eq. (4.6). Notice 
that the last three terms in Eq. (4.7) will cancel out. On the other hand in the tight coupling limit expansion one 
can set v^ 1 — v^ 1 and v^ 1 — v^ 1 in the remaining terms on the R.H.S. of Eq. (4.7). Thus Eq. (4.6) becomes 



l + i2 uu 1 + i? v ' 3(1 + J2) W " M ' J 7 1 + i? 

' -($« + tfWjA^ - -tL-Pv™ 2 - 2-^-($« + - /JL^ 1 ) („(i>* - , (4.8) 



2(l + i?) v y uu 1 + 7 l + i? v y 1 + i? 

where in the tightly coupled limit we are neglecting the first-order quadrupolc and (higher-order moments) of the 

photon distribution since it is suprressed by 1/r with respect to the other terms. Next for the term like t i 8^\v^ 1 — 

) we employ the relation previously derived in Eq. (4.3) with 5^ — 3Aqq*/4 and we use the first order tight 
coupling equations (3.1) and (3.8) in order to further simplify Eq. (4.8). We finally obtain 

R 1 A (2) < 4 

(2)i/ , n n (2)» , i^ao_ , $(2),i = S i (4 Q) 

7 +rt l + i?, 7 + 41 + i? + * v ' 1 yj 

where 

6y " 4(1 + 4 "l + i? +2 (l + i?)2 A ° ^ + 4(l + i?)2 l A °° J + 3(l + i?) 7 dj 7 

+ 'ifi?^ 1 ^^ - 2(<i>(1) + * (1))<i>(1) ' 1 - 2(ttr) (<&(1) + * (1) ) A o 1 o ),< + 3(TTe) (?? - + 

n fy v (m |-A-A« f -^-^ - ]^-) , (4.10) 



1 + i? 7 21 + i? 00 ^l + i? 7 41 + i? / 

We have spent some time in giving the details of the computation for the photon Boltzmann equations at second-order 
in the perturbations. As a summary of the results obtained so far we refer the reader to Eqs. (3.1)-(3.8) and Eqs. (4.4)- 
(4.9) as our master equations which we will solve in the next sections. In particular Eq. (4.9) is the second-order 
counterpart of Eq. (3.8) for the photon velocity in the tight coupling regime. Notice that there are two important 
differences with respect to the linear case. One is that, in Eq. (4.9), there will be a contribution not only from scalar 
perturbations but also from vector modes which, at second-order, are inevitably generated as non-linear combinations 
of first-order scalar perturbations. In particular we have included the vector metric perturbations uj % in the source 
term. Second, and most important, we have also kept in the source term the second-order quadrupolc of the photon 
distribution iT^ u . At linear order we can neglect it together with higher order moments of the photons since they 
turn out to be suppressed with respect to the first two moments in the tight coupling limit by increasing powers 
of 1/r. However in the next section we will show that at second order this does not hold anymore, as the photon 
quadrupole is no longer suppressed. 

Finall following the same steps that lead to Eq. (3.9) at linear order we can derive a similar equation for the 
second-order photon energy density perturbation Aqq which now will be characterized by the source terms Sa and 

si, 

(a$" - 4*( 2 )") + (a$' - 4*( 2 )') - c 2 v 2 (a&> - 4*( 2 )) = iv 2 (V 2 ) + ^) + s' A + n^s A 

- piS l v. (4.11) 



C. Second-order quadrupole moment of the photons in the tight coupling limit 

Let us now consider the quadrupole moment of the photon distribution defined in Eq. (3.4) and show that at 
second-order it cannot be neglected in the tightly coupled limit, unlike for the linear case. We first integrate the 
R.H.S. of Eq. (2.2) over dQ n (n l n-? — £ u /3)/47r and then we set it to be vanishing in the limit of tight coupling. 
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The integration involves various pieces to compute. For clarity we will consider each of them separately. The term 

(2) 
00 



[2] 

A 00 does not contribute. For the third term we can write, from Eq. (2.4) 



1 ^ 2 V4^ Af 2) A (2) 1 ^ .(2), sf I 4tt 

- 2 E h^ A ^ Y2m = ~W ~ ToE E A ®H) V^TT r£m ' (4 - 12) 

m=-2 t±1m=-t 

so that the integral just brings IT^^/IO, since the only contribution in Eq. (4.12) comes from A^ 2 ' /10 with all the 
other terms vanishing. The following nontrivial integral is 



/ 



*J (W _ I^A A «P 2 (v • n) = fi fc fi, / ^ (W - |*«) A« QnV - ±) = (W - ^ ) (4.13, 



where the baryon velocity appearing in P2 (v • n) is first order and we make use of the following relations 

J dn n l = J dtt n l n ] n k = , J ^ nW = ^S ij , J ^ nVnV = ^(S ij 6 kl + S ik S lj + S il S jk ) . (4.14) 
The integrals of S^A^ , #(vn) and v< 2 ) • n vanish and 

vA^ J ^ 2 (v • n) (nV - = ^vA^ - = ^ (vV - \s ij v 2 ^j , (4.15) 

where in the last step we take A^ = 4w/3 in the tight coupling limit. Similarly the integral of 14(v • n) 2 brings 

lW J ^n k n e (nW - = ^ (v i v j - \^v 2 ^j . (4.16) 

The integral of 2(v • n)AW can be performed by expanding the linear anisotropics as A^ = J2 l (2£+ l)A^P^(v • n). 
We thus find 

2v h %m J f^ n3 _ Iji^ nmA (D + Q{1 > 2) = 16 ^ v i v j _ l&v 2 ^ , (4.17) 

where we have used Eq. (4.14) and 0(£ > 2) indicates all the integrals coming from the multipoles £ > 2 in the 
expansion (for I — and 1 = 1 they vanish.) In fact we have dropped the 0(£ > 2) since they are proportional 
to firts-order photon moments £ > 2 which turn out to be suprressed in the tight coupling limit. Finally the term 

proportional to (v • n)A 2 ^(l — P2(v • n)/5 gives a vanishing contribution. 

Collecting all the various pieces we find that the third moment of the R.H.S. of Eq. (2.2) is given by 



(4.18) 



Therefore in the limit of tight coupling, when the interaction rate is very high, the second-order quadrupole moment 
is given by 

n (2)y ~ § - , (4.19) 

by setting Eq.(4.18) to be vanishing (the term multiplying si^ goes to zero in the tight coupling limit since it just 
comes form the first-order collision term). At linear order one would simply get the term QrTI^^/lO implying that, 
in the limit of a high scattering rate r', 11^ ^ goes to zero. However at second-order the quadrupole is not suppressed 
in the tight coupling limit becasue it turns out to be sourced by the linear velocity squared. 



V. SECOND-ORDER CMB ANISOTROPIES GENERATED AT RECOMBINATION 



The previous equations allow us to follow the evolution of the monopole and dipole of CMB photons at recombi- 
nation. As at linear order, they will appear in the expression for the CMB anisotropies today A^ 2 ^ (k, n, 770) together 
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with various integrated effects. Our focus now will be to obtain an expression for the second-order CMB anisotropies 
today A^ 2 ^ (k, n, 770) from which we can extract all those contributions generated specifically at recombination due to 

(2) 

the non-linear dynamics of the photon-baryon fluid. This expression will not only relate the moments A)^ today to 
the second-order monopole and dipole at recombination as it happens at linear order, but one has to properly account 
also for additional first-order squared contributions. Let us see how to achieve this goal in some details. 

As in linear theory (see e.g. [33, 34]) it is possible to write down an integral solution of the photon Boltzmann 
equation (2.2) in Fourier space. Following the standard procedure for linear perturbations, we write 



A®' + ikfxA^ - r'A^ = e -^+ra ( 2)eikm = ^ 

dr\ 



in order to derive a solution of the form 

A( 2 )(k,n,r ?0 ) 



no 



drjS(k 7 n, r])e 



ikn(r]-T]o) 



Here \i = cos-d = k • n is the polar angle of the photon momentum in a coordinate system such that e3 
second-order the source term has been computed in Ref. [22] and can be read off Eq. (2.2) and Eq. (2.6) to be 

S = -t'A$ - An^f + 4*( 2 >' + 8A«(4-«' - n^f) - 2n 4 ($« + *«)(A« + 4$«), i 
r i i^AW 



(5.1) 
(5.2) 

k. At 



— T 



-\ t ^A^r 2m („) + 2^)(A«-AW + 4v.n + iAWp 2 (v.n)) + 4v( 2 ).n 



+2(v • n) 



A 



(i) 



3A 



(i) 



A. 



(i) 



1 



-P 2 (v-n) 



- vA{ 1} (4 + 2P 2 (v • n)) + 14(v • n) 2 



2v 2 



(5.3) 



The key point here is to isolate all those terms that multiply the differential optical depth r'. The reason is that in 
this case in the integral (5.2) one recognizes the visibility function g(rj) = —e~ T r' which is sharply peaked at the time 
of recombination and whose integral over time is normalized to unity. Thus for these terms the integral just reduces 
to the remaining integrand (apart from the visibility function) evaluated at recombination. The standard example 

that one encounters also at linear order is given by the first term appearing in the source S, Eq. (5.3), that is — t'Aq . 
The contribution of this term to the integral (5.2) just reduces to 

rVo 

A( 2 )(k,n,r /0 ) = / d V e lk ^>-^e- T (-T')A$ ~ e^^-^A^*) , (5.4) 
Jo 

where r]* is the epoch of recombination and, in the multipole decomposition (2.5), Eq. (5.4) brings the standard result 

A£>( % ) <x *$(V.)MKV. - %)) , (5.5) 

having used the Legendre expansion e lk x = ^2 £ (iY(2£+l)jg(kx)P£(k-x). In Eq. (5.5) the monopole at recombination 
is found by solving the Boltzmann equations Eqs. (4.4)-(4.9) derived in tight coupling limit. 

Looking at Eq. (5.3) we recognize immediately some terms which multiply explicitly r' (the first one discussed in 
the example above and the last two lines of Eq. (5.3)). However it is easy to realize from the standard procedure 
adopted at the linear-order that such terms are not the only ones. This is clear by focusing, as an example, on the 

(2) 

term — 4n J <i> \ in the source S which appears in the same form also at linear order. In Fourier space one can replace 
the angle ^ with a time derivative and thus this term gives rise to [33, 34] 

) e ->$( 2 ) = -4/ dri^e-^t 
Jo d V v 



- Aik 



dr\ e 



ikn(v-Vo). 



Jkfj.(r]-r]o) 



dr)e ik ^-vo) e -T ^(2)'_ T / $ (2; 



(5.6) 

where, in the last step, we have integrated by parts. In Eq. (5.6) the time derivative of the gravitational potential 
contributes to the Integrated Sachs- Wolfe effect, but also also a r' results implying that we have also to evaluate $^ 2 - ) 
at recombination. Thus, in the following, we look for those terms in the source (5.3) which give rise to a r' factor in 

the same way as for — 4n*$ \ . In particular let us consider the combination in Eq. (5.3) 

C = SAW^W - n l ^ ] ) - 2n l (^ + *«)(AW + = SA^W - Sr^AW 1 )),; + 4$( 1 WA^ 1) 

- 8n 4 ($(i)2) ji7 (5.7) 
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where for simplicity we are setting ~ We already recognize terms of the form n l di{-). Moreover we can use 

the Boltzmann equation (2.1) to replace n l A { l ] in Eq. (5.7). This brings 



C = 8AW 1 )' - 4*«A< 1 ) / - 16^^' - 8n 4 (A«$«), - 16n J ($« 2 ) , t 
4t'*« [a$ - A« + 4v« • n + iA 2 1} P 2 (v • n) 



(5.8) 



In fact we will not be interested for our purposes in the first three terms of Eq. (5.8), since they will not contribute 
to the anisotropies generated at recombination. 

Therefore, as a result of Eqs. (5.3), (5.6) and (5.8), we can rewrite the source term (5.3) as 



S = S* + S' 



(5.9) 



where 
S* - 



A$ + ^ ~ \ E + 2«W (A« - A« + 4v • „ + \ A«P 2 (v . „)) + ^ ■ n 



+2(v • n) 



AW+SA^-AW (l-5p 2 (v.n) 



- «Ai 1} (4 + 2P 2 (v • n)) + 14(v • n) 2 - 2v 2 + 8A (1) $ (1) 



16$ (1)2 + 4^ (1) 



A« - A« + 4v« • n + ^A 2 1} P 2 (v • n) 



(5.10) 



and 



S' = 4($( 2 ) + f( 2 ))' - 8^n J - 4x- j n i n j - 2 [($« + ^'^jTiW - ($« + tfW)- 4 

+ 8(aW$W)' + saW^w - 4fWAW + le*^)*^)' . 



gAW 

dri 1 



(5.11) 



In Eq. (5.9) 5* contains the contribution to the second-order CMB anisotropies created on the last scattering surface 
at recombination, while S' includes all those effects which are integrated in time from the last scattering surface up 
to now, including the second-order Integrated Sachs- Wolfe effect and the second-order lensing effect. Since the main 
concern of this paper is the CMB anisotropies generated at last scattering, from now on we will focus only on the 
contribution from the last scattering surface S* . 



A. Multipole moment decomposition 



( 2") 

The expression of the photon moments A^ can be obtained from Eq. (2.5). Such a decomposition can be achieved 
by first expanding the source term S as 



e m= 



4ir 
2i+l 



Y em (n) 



(5.12) 



and then taking into account the additional angular dependence in the exponential of Eq. (5.2) by recalling that 

e ,k.x = J2( t ) e (2e + l)j t (kx)P t (k ■ x) . (5.13) 

e 

Thus the angular integral (2.5) just reduces to compute the expansion coefficients of the source term 

t-2 



a^(Mo) = {-i)- m {-i)- l (2e + i) / d V e-^J2 E H)' 2 ^E^ Wi + mAKv-m)) 

J ° £ 2 m 2 = -l 2 



I \ m 2 -m 



(5.14) 
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where the Wigner 3 — j symbols appear because of the Gaunt integrals 

GuZ 2 " 13 = J d 2 nY limi (h)Y l2m2 (n)Y h m 3 (A) 

/ (2Z 1 + 1)(2Z 2 + 1)(2Z 3 + 1) (lx t 2 4 \( ti h is \ 
V 4tt V J \ mi m 2 m 3 J 

The observed anisotropies generated at the last scattering surface come from the source term containing a — r' 
factor: this allows to solve the time integral in Eq. (5.14) by evaluating the integrand at rj = ry* given that the visibility 
function g(rj) = — T'e~ r is peaked at the time of recombination. 



VI. TIGHTLY COUPLED SOLUTIONS FOR THE SECOND-ORDER PERTURBATIONS 

In this section we will solve the tightlty coupled limit of the Boltzmann equations (4.4) and (4.6) at second-order in 
perturbation theory. We will proceed as for the linear case, focusing on the two limiting cases of perturbation modes 
entering the horizon respectively much before and much after the time of equality. The solution of Eq. (4.11) can be 
written as 

[1 + i?(r 7 )] 1 / 4 (A( 2 ) - 4*( 2 )) = Acos[kr s (ri)} + B sin[fcr s (77)] 

- 4 A £ drf{1 + RtfyfH (>)(,/) + sin[fc(r a („) r.tf))] 

+ X [ dV ' [l + i?(?? ' )]3/4 { S ' A + TTR SA t ikiS v) si ^ k (r s (v) r s ( V '))} ,(6.1) 

where the source terms are given in Eq. (4.5) and (4.10). Notice that we can write S' A + j^Sa — (S^(l + R))' /1 + R 
so that we can perform an integration by parts in Eq. (6.1) leading to 

AT 

[1 + #(??)] ^(Aoo -4* (2) ) = ,4cos[fcr s (77)] +B sin[fcr s (?7)] - ^-5 A (0) sm[kr s (r,)} 

4-|^'^ [1 + W] 3/4 (^(2) {v , ) + ^m^ si n[k(r s (r,)-r s ( V '))} 

+ r H S A ( V ') (1 + R( V ')) 1/4 co S [k(r s (r,) - r a (rf))] 
Jo 

l^jf" ( | J| '5Mn')(l + %')) 3/4 S m[k(rM-rs(v'))} 
+ ^ J\ V 'SA(i){l + R(ri'))- l/i R\^)MKrs^)-r s (i))]. (6.2) 



A. Setting the initial conditions: primordial non-Gaussianity 

The integration constants A and B are fixed according to the initial conditons for the second-order cosmological 
perturbations. These refer to the values of the perturbations on supcrhorizon scales deep in the radiation dominated 
period. We will consider the case of initial adiabatic perturbations, for which there exist some useful conserved 
quantities on large scales which as such carry directly the information about the initial conditons. 

In the standard single-field inflationary model, the first seeds of density fluctuations are generated on super-horizon 
scales from the fluctuations of a scalar field, the inflaton [1]. Recently many other scenarios have been proposed 
as alternative mechanisms to generate such primordial seeds. These include, for example, the curvaton [35] and 
the inhomogeneous reheating scenarios [10], where essentially the first density fluctuations are produced through 
the fluctuations of a scalar field other than the inflaton. In order to follow the evolution on super-horizon scales of 
the density fluctuations coming from the various mechanisms, we use the curvature perturbation of uniform density 
hypersurfaces ( = + /2 + •••, where = — "J'W — HSp^ jp' and the expression for is given by [36] 

C (2) = _^(2) _ H S _^_P + AC (2) f (6 3) 
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with 



AC ( 2 ) = 2H 5 —P- 6 —P + 2 6 —P(¥^> + 2U¥^) ( S —f-) (h P --W- 2<yA + 2*« 2 . 

p' p' p' \ p J \ p J 



(6.4) 



The crucial point is that the gauge-invariant curvature perturbation £ remains constant on super-horizon scales after it 
has been generated during a primordial epoch and possible isocurvature perturbations are no longer present. Therefore, 
we may set the initial conditions at the time when £ becomes constant. In particular, provides the necessary 
information about the "primordial" level of non-Gaussianity generated either during inflation, as in the standard 
scenario, or immediately after it, as in the curvaton scenario. Different scenarios are characterized by different values 

of (( 2 \ For example, in the standard single-field inflationary model = 2 (C*- 1 ^) + O (e, 77) [5, 37], where e and 
r] are the standard slow- roll parameters [1]. In general, we may parametrize the primordial non-Gaussianity level in 
terms of the conserved curvature perturbation as in Ref. [38] 



C (2) =2a NL (cW) 2 , (6.5) 



where the parameter cinl depends on the physics of a given scenario. For example in the standard scenario onl — 1, 
while in the curvaton case onl = (3/4r) — r/2, where r ~ {pa/ p)t> is the relative curvaton contribution to the total 
energy density at curvaton decay [4, 8]. In the minimal picture for the inhomogeneous reheating scenario, onl = 1/4. 
For other scenarios we refer the reader to Ref. [4]. One of the best tools to detect or constrain the primordial large- 
scale non-Gaussianity is through the analysis of the CMB anisotropies, for example by studying the bispectrum [4]. 
In that case the standard procedure is to introduce the non-linearity parameter /nl characterizing non-Gaussianity 
in the large-scale temperature anisotropies [4, 26, 30]. To give the feeling of the resulting size of /nl when |onl| S> 1, 
/nl 5a NL /3 (see Refs. [4, 38]). 

The conserved value of the curvature perturbation Q allows to set the initial conditions for the metric and matter 
perturbations accounting for the primordial contributions. At linear order during the radiation-dominated epoch and 
on large scales — — 2\I>( 1 )/3. On the other hand, after some calculations, one can easily compute A(( 2 ) for a 
radiation dominated epoch 

AC (2) ^(* (1) ) 2 , (6.6) 

where in Eq. (6.4) one uses that on large scales 5^p^/p y = —2^^ and the energy continuity equation 5^' p 1 + 
mS^p 7 - 4*( 1 )' ( o 7 = 0. Therefore we find 

C (2) = _^(2) + ^ + 7^(l) 2(0)5 (6?) 

where we are evaluating the quantities in the large scale limit for 77 — > 0. Using the parametrization (6.5) at the initial 
times the quantity — 4 , 3>( 2 ) is given by 

- 4* (2) = 2(9a NL - 7)* (1)2 (0) . (6.8) 

Since for adiabatic perturbations such a quantity is conserved on superhorizon scales, it follows that the constant 
B = and A = 2(9a NL - 7)¥^ 2 (0). 

Eqs. (6.1) and (6.2) are analytical expressions describing the acoustic oscillations of the photon-baryon fluid induced 
at second-order for perturbation modes within the horizon at recombination. In the following we will adopt similar 
simplifications already used for the linear case in order to provide some analytical solutions. In particular, if in 
Eq. (6.2) we treat R as a constant we can write, using the initial conditions determined above, 

(A$ - 4*< 2 )) = 2(9a N L - 7)*< 1 ) 2 (0) cos[fcr s (77)] - ^yS A (0) sm[kr s ( V )} 



dr f ($( 2 )(r/) + *( 2 )(7/)) sm[k(r s (n) - r s (v'))} 



4 k 
3 c s 

+ ^ d v 'S A {n') cos[fc(r 8 (»7) -r.fa'))] 
Jo 



Vlt s [ H ^W*?)-^'))]. (6.9) 

Notice that we have also dropped the occurence of R in $< 2 ) +*< 2 )/i?. 
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VII. PERTURBATION MODES WITH k > k, 



In order to study the contribution to the second-order CMB anisotropics coming from perturbation modes that 
enter the horizon during the radiation dominated epoch, we will assume that the second-order gravitational potentials 
are the ones of a pure radiation dominated universe throughout the evolution. Though not strictly correct, this 
approximation will give us the basic picture of the acoustic oscillations for the baryon-photon fluid occuring for these 
modes. Also for the second-order case, in Section IX we will provide the appropriate corrections accounting for the 
transition from radiation to matter domination which is indeed (almost) achieved by the recombination epoch. Before 
moving into the details a note of caution is in order here. At second order in the perturbations all the relevant 
quantities are expressed as convolutions of linear perturbations, bringing to a mode-mode mixing. In some cases 
in our treatment for a given regime under analysis (k » k eq or k <C k eq ) we use for the first-order perturbations 
the solutions corresponding to that particular regime, while the mode-mode mixing would require to consider in 
the convolutions (where one is integrating over all the wavenumbers) a more general expression for the first-order 
perturbations (which analytically does not exist anyway). For the computation of the CMB bispectrum this would be 
equivalent to consider just some specific scales, i.e. all the three scales involved in the bispectrum should correspond 
approximately to wavenumbers k 3> k eq or k <C k eq , and not a combination of the two regimes (a step towards the 
evaluation of the three-point correlation function has been taken on Ref. [27] where it was computed in the in so-called 
squeezed triangle limit, when one mode has a wavelength much larger than the other two and is outside the horizon). 

Having learned that at linear order the regime k 3> k eq can be solved in the alternative way described by Eq. (3.27) 
and (3.29), we adopt the same procedure can be adopted also at second-order: we will use Eq. (4.4) where we can 
neglect the gravitational potential term , 3>( 2 )'. The reason is again that also the second-order gravitational potentials 
decay at late times as t/~ 2 , while the second-order velocity oscillate in time. Let us now see that is some details. 

The evolution equation for the gravitational potential 'J'' 2 ) is given by Eq. (B.18) and is characterized by the source 
term S 1 , Eq. (B.19). In particular the source term contains the second-order quadrupole moment of the photons 

4 2)y . We saw in Section IV C that at second-order the quadrupole moment is not suppressed in the tight coupling 
limit, being fed by the non-linear combination of the first-order velocities, Eq. (4.19). For the pertubation modes we 
are considering here the velocity at late times is oscillating being given by Eq. (3.27) in Fourier space. Since the linear 
gravitational potential (3.26) decays in time and for a radiation dominated period Ti = 1/r], it is easy to check that 
the dominant contribution at late times to the source term S* 7 simply reduces to 

5 7 * ^^nf , = F ^ ^ C <> (0)<> (0) sin( fclCs?7) singed , (7.1) 



where we have used Eq. (3.27), 



cjkik 2 



and the sound speed is c s = 1/ \/3(l + R). Before proceeding further let us explain the notation that we are using. 
The equivalence symbol will be used to indicate that we are evaluating the expression in Fourier space. At second- 
order in perturbation theory most of the Fourier transforms reduce to some convolutions. We will not indicate these 
convolutions explicitly but just through their kernel. For example in Eq. (7.1) by F(ki,k 2 ,k) we actually indicate 
the convolution operator 

y*d 3 M 3 fc2<5 (1) (ki +k 2 -k)F(k 1; k 2 ,k). (7.3) 
In the specific case of Eq. (7.1) the kernel is given by 

f (k 1 ,k a ,k) = (k - k ^ k - k2) -lk 1 .k 2 . (7.4) 

The choice of these conventions is due not only for simplicity and to keep our expressions shorter, but also because at 
the end we will be interested to the bispectrum of the CMB anisotropics generated at recombination, and the relevant 
expressions entering in the bispectrum are just the kernels of the convolution integrals. 

Having determined the leading contribution to the source term at late times, we can now solve the evolution 
equation (B.18). Since the source term scales like r\~ 2 it is useful to introduce the rescaled variable 
Eq. (B.18) then reads 



x "+[k 2 cl-^) X = V 2 S 1 . (7.5) 
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For perturbation modes which are subhorizon with k-q » 1 the solution of the homogeneous equation is given by 

Xhom. = Acos(kc s j]) + B sin(c s kr]) , (7.6) 
from which we can build the general solution 



X Xhom. "T X+ 







(7.7) 



where = — kc s is the Wronskian, and \+ = cos(fcc s r;), x_ = sin(c s fc?7). Using Eq. (7.1) the integrals involve 
products of sines and cosines which can be performed giving 

FC m T m k [2fc 1 fc 2 cos(fciC s r?)cos(fc 2 c s r?) - 2k 1 k 2 cos(kc s r)) + (k\ + k\ - k 2 ) sin(fc 1 c s ?7) sin(fc 2 c s ry)l 

X = Xhom.-- 2 -*k 1 (°)*k 2 (°)— 



7-4 _i_ 7,4 _j_ 7,4 9 £,2 £.2 97,2 7,2 0h 2 h 2 



Thus the gravitational potential \I/( 2 ) at late times is given by 



(7.8) 



*k kn) = -3* (2) (o) 



cos(kc s rj) 



FC 

V 2 c 2 s 



{kc s -q) 2 



[2fcifc 2 cos(fcic ;j ?7)cos(fc2C s ?7) — 2fcifc 2 cos(fcc s 77) + (kf + k 2 — k 2 ) sin(kic s r]) sin(k2C s r])] 



U4 I 1,4 I 1,4 _ OjU2l,2 _ 97,2 7,2 _ 97,27,2 



(7.9) 



where we have set the integration constant _B = and A — — 3\E , ^ 2 - ) (0)/(fcc s ) 2 in order to match the homogenoeus 
solution at late times which has the same form as Eq. (3.26). Here <J>( 2 )(0) is the intial condition for the gravitational 
potential taken on large scales deep in the radiation dominated era which will be determined in Section VII B. 

Eq. (7.9) shows the result that we anticipated: also at second order the gravitational potential varies in time 
oscillating with an amplitude that decays as r\~ 2 . Let us then take the divergence of the (i — 0) Einstein equation (A. 10) 
expanded at second-order 



di 



2H 2 d l 



2 V i 



P)< + ($W + fW) B f 



+ A (1 V 1}i 



which, using the first-order (i — 0) Einstein equation (B.17) and ~ "J'W, reduces to 



2 2 



-2H 2 d l 



2 V l + ^00 U 7 



(7.10) 



(7.11) 



Since during a radiation dominated period is given by Eq. (3.26) at late times, it is easy to see that (vI/W'c^'I'M) 
will be oscillating and decaying as r]~ 4 and thus can be neglected with respect to v[/( 2 )' ; which oscillates with an 
amplitude decaying as r)~ 2 . Also Tt®^ turns out to be subdominant. Recall that $^ 2 ' = *( 2) - (sec Eq. (A. 12)) 
and is dominated by the second-order quadrupole of the photons in Eq. (B.19), so that <J>( 2 ) scales like *^ 2 ^ but 
there is the additional damping factor of the Hubble rate H = l/r]. Thus the dominant terms give 



.^v^'^a^W) 



(7.12) 



Eq. (7.12) is the equivalent of Eq. (3.28) and it allows to proceed further in a similar way as for the linear case by 
using the results found so far, Eqs. (7.9) and (7.12), in the energy continuity equation (4.4). In Eq. (4.4) the first- and 
second-order gravitational potentials can be neglected with respect to the remaining terms given by Aqq and v^ 1 
which oscillate in time. Thus, replacing the divergence of the second-order velocity by the expression (7.12), Eq. (4.4) 
becomes 



A (2)/ - 



_v 2 *( 2 )' + - 
m 2 3 



dv^A (1) + fA (1)2 V 



(7.13) 



which, using the first-order equation (3.1), further simplifies to 



v (2)' 
*00 



-v 2 *< 2 )' 

1 5 



(7.14) 
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where we have kept only the dominant terms at late times. 

The gravitational potential is given in Eq. (7.9), so the integration of Eq. (7.14) gives 

A$ = 6* (2) (0)cos(/cc s r?) 

FC (i) (i) 2 [2fcifc 2 cos(fcic s 77)cos(fc 2 c s 77) - 2k 1 k 2 co$,(kc s 'q) + (k\ + fc| - k 2 ) sin{k 1 c s 'q) sin(fc 2 c s ?7)] 

UO^ ft ^ I A' Zj ft ft 2 ZjAj^ ft* Zift^A' 

(7.15) 

Needless to say, modes for k 3> fcr>, where fc^ 1 indicates the usual damping length, are supposed to be multiplied 
by an exponential e~C c / fcl >) (see, e.g. [33]). 



A. Vector perturbations 

So far we have discussed only scalar perturbations. However at second-order in perturbation theory an unavoidable 
prediction is that also vector (and tensor) perturbation modes are produced dynamically as non-linear combination of 
first-order scalar perturbations. In particular notice that the second-order velocity appearing in Eq. (5.10), giving rise 
to a second-order Dopplcr effect at last scattering, will contain a scalar and a vector (divergence free) part. Eq. (7.12) 
provides the scalar component of the second-order velocity. We now derive an expression for the velocity that includes 
also the vector contribution. 

The (second-order) vector metric perturbation oj 1 when radiation dominates can be obtained from Eq. (B.22) 

- \ W + 3«V = -AH 2 (a* - *|) ^ + A«<)^ , (7.16) 

where we have dropped the gravitational potentials \f'( 1 ) ~ which are subdominant at late times. On the other 
hand from the velocity continuity equation (4.9) we get 

v^' + \ = \ Ko )2 ) + \^%v^ U» \d k ^ , (7.17) 

neglecting the term proportional to R and the decaying gravitational potentials. Using the tight coupling equations 
at first-order (3.1) and (3.2), and integrating over time one finds 

„(2)i + 2(vWA$) = -2^ \ J drf*®* - \ J d v 'd k n^ . (7.18) 

We can thus plug Eq. (7.18) into Eq. (7.16) to find that at late times (for ki] > 1) 

VV = -m 2 (S) - J d V 'd k U^ . (7.19) 

We will come later to the explicit expression for the term on the R.H.S. of Eq. (7.19). Here it is enough to notice 
that the second-order quadrupole oscillate in time and thus u l will decay in time as H 2 = l/rj 2 . This shows that uo l 
in Eq. (7.18) can be in fact neglected with respect to the other terms giving 



„(2)i 



-2(,WA«) - \ J drj'A^ \ J d v 'd k U^ . (7.20) 



It can be useful to compute the combination on the R.H.S. of Eq. (7.19) (5j — d l dj /V 2 )d k TL^ kj . The second-order 
quadrupole moment of the photons in the tightly coupled limit is given by Eq. (4.19), and 

dkHVW = | [d k (v k v j ) - 2v k &v k ] = I [v j 8 k v k - v k &v k ] , (7.21) 
where in the last step we have used that the linear velocity is the gradient of a scalar perturbation. We thus find 

U 3 - d k U^ = | (v*d k v k - v k d l v k ) - ||1 [(8 k v k ) 2 + v^d 3 d k v k - d 3 v k &v k - v k V 2 v k ] . (7.22) 
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Notice that if we split the quadrupole moment into a scalar, vector (divergence-free) and tensor (divergence-free and 
traceless) parts as 

n (2)kj = jjpj.fcj _ 1 V 2 6 kj U (2) + n (2)k,j + n (2),,fe + n PW ; (7 23) 

o 



5% - ^ ) 8 k U^ = V 2 nf J , (7.24) 



then it turns out that 

where H^ 1 is the vector part of the quadrupole moment. Therefore one can rewrite Eq. (7.19) as 

lo 1 = -3H 2 J d v 'n^ . (7.25) 

B. Initial conditions for the second-order gravitational potentials 

In order to complete the study of the CMB anisotropics at second-order for modes k » k eq we have to specify the 
initial conditions ^^(O) appearing in Eq. (7.15). These are set on super-horizon scales deep in the standard radiation 
dominated epoch (for 77 — > 0) by exploiting the conservation in time of the curvature perturbation £. On superhorizon 
scales C^ 2 ^ is given by Eq. (6.7) during the radiation dominated epoch and, using the (0 — 0)-Einstein equation in the 
large scale limit = -2$ (2) + 4<f> (1 ) 2 , we find 

C (2) = -^*( 2 )(0) - i (V 2 )(0) - *^(0)) + ^* (1)2 (0) . (7.26) 

The conserved value of is parametrized by = 2a^(^ 2 , where, as explained in Section VI the parameter 
&nl specifies the level of primordial non-Gaussianity depending on the particular scenario for the generation of 
the cosmological perturbations. On the other hand at second-order the gravitational potentials differ according to 
Eq. (A. 12), which for superhorizon modes during radiation domination gives 

$( 2 )(0)-* (2) (0) = -Q (2) (0), (7.27) 

where 

q( 2 >(o) = -2v- 2 a fc $(i)(o)a fc $( 1 )(o) + 6^ (^$«(o)d,$«(o)) + ^ 2 ^nf; , (7.28) 

where we are evaluating Eq. (B.20) in the limit kr\ <C 1. The gravitational potential (B.16) just reduces to the constant 
$^^(0), while the contribution from the second-order quadrupole moment in this limit reads 

^ 2 ^<); = Jjjf** („V - f = -3|^< 1 i )(0)< 1 2 ) (0)sin(. lCs ,)sin( fc2Cs ,) - §F^(0)^(0) , 

(7.29) 

where F and C are defined in Eqs. (7.3) and (7.2). Therefore we find that in Fourier space 

Q (2) (0) = 33^ 1 i ) (0)< ) (0), (7.30) 
and from Eq. (7.26) we read off the intial condition as (convolution products are understood) 

*%m£(o)- (7-3i) 



*(2)(0) = 



- 3 (a NL -l) + ll F < k " k - k > 



k 2 

C. Multipole moments 



In this Section we give the expression for the CMB multipole moments observed today which are due to the 
perturbations of the photons at the last scattering surface. Therefore we make use of Eq. (5.14) where we just 
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consider the part 5* of the source term. As explained in Sec. V S* contains a — r' factor which reduces the time 
integral in Eq. (5.14) by evaluating the integrand at 17 = 77* given that the visibility function g(jf) = — t' e~ T is peaked 
at the time of recombination. Therefore we evaluate 5* at recombination in the tightly coupled limit for the modes 
k ^> k eq using the previous results and we decompose it according to Eq. (5.12). 

First we use the solution for the photon-baryon fluid of Eq. (7.20) in Eq. (5.10) to find 



A oo } - \ E ^A^y 2m( n) - n • V J */A# - 3n • vj drfU™ 2(v • n)A« + 2(v • n)A (1 



m=-2 

- uA^(4 + 2P 2 (v • n)) + 14(v n) 2 - 2u 2 



(7.32) 



Notice that in Eq. (5.10) we have neglected all the terms depending on the gravitational potentials since they decay 
in time, the terms proportional to the linear dipolc which is suppressed in the tight coupling limit, and the terms 
proportional to (A ^ — A' 1 ) + 4v • n) which is suppressed being just the first-order collision term. 

For the decomposition of S* in multipolc moments Adjust gives A$ 5eo5 m o. Similar terms, which do not carry 
angular dependence in Eq. (5.10), are — 2v 2 and — 4vA^\ Notice that in the limit of tight coupling we can use 
A^ = 4w/3. For the terms which arc quadratic in the velocities it is convenient to write 



14(v • n) 2 - 2v 2 - ivA^ = \^{n l v? - T S lj )v t v 



(7.33) 



The term 14(n l n : ' — ^8 l: >)v i Vj can be decomposed with multipoles given by (in Fourier space) 



14H)" 



21- 



47T 



1 /4tt 
3" 



w(ki>(k 2 )(-i)- m j2 n* mi (ki)r* lm2 (k 2 )^r 



1 



:V 2 S m S m0 , 



(7.34) 



where in Fourier space, for the first-order velocity we use the convention 

v(ki) = iw(ki)ki , 



(7.35) 



and the convolution products are implicitly assumed in a similar way as in Eq. (7.3). In order to derive Eq. (7.34) 
and the following expressions we use the addition theorem of the spherical harmonics 



(7.36) 



Notice that the term — 4wA^ 1 - ) P 2 (v • n) can be written as — 4(r 



S lJ )viVj 



of the same type as that in Eq. (7.34). 



The multipoles of (-2A{,J ) ) v • n) are 2i Agg'' (k 2 )u(ki)y l ^ Tl (ki)(5i£ using the same rules as above. 
Let us now consider the term 2(v • n)A( 1 ^. From Eq. (2.4) we can write 



(7.37) 



where we are neglecting higher-order multipoles in the tight coupling limit. Therefore the multipoles of 2(v • n)A^ 1 ) 
are given by 



2i(-i) 



2£+l 

47T 



47T 



V (k!)- ]T v* mi (U) 



A$ (k 2 )6 n S mmi + AY' (k 2 )3W ^(-1) 



(!)/ 



(7.38) 



Notice that for £ = and m = Eq. (7.38) gives 8v 2 /3 which then will cancel the second term on the R.H.S. of 
Eq. (7.34). This can be accounted for by simply neglecting such a term in Eq. (7.34) and writing Eq. (7.38) by 
specifying f ^0,m^0. 

The term — n • V f drj'A^ has expansion coefficient 



(7.39) 
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Finally the expansion coefficients for the second term in Eq. (7.32), - Y^m=-2 v / 47rA^y 2 m(n)/(25 3/2 ) 
to while the term -3n • V / dr)'ltf )ij has expansion coefficients 



-k 2 + ki • k 2 



where we have used Eq. (4.19). 

Collecting all the previous results we find 



~>*lm 



-A (-i)- e -' 2i+1 



Air 



4nA$6 m 5 m0 + J 8 7 r A /^ U (k 1 )Ai^(k 2 )(-l)- m g™f ~ m (£ + 0;m ? 0) 



d? 7 ' w (k 1 )«(k 2 )n(k 1 )(5, 1 , 



reduces 



(7.40) 



io W (k 1 >(k 2 )(-i)- m £ K 1 * roi (k 1 )y 1 * ma (k 2 )g 



mini2—m 
lit 



+ S n 5 mQ k I dri'A$(Ti') 



-k 2 + ki • k 2 



drj'v(ki)v(k2)Y 1 * m (ki)S£i + S e2 : 



i (2) 

*2m 
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(7.41) 



Eq. (7.41) is all we need to get the multipole moments today given by Eq. (5.14). 



VIII. PERTURBATION MODES WITH k < k e 



Let us consider the photon perturbations which enter the horizon between the equality epoch and the recombination 
epoch, with wavelenghts rj^ 1 < k < T)~g. In fact, in order to find some analytical solutions, we will assume that by the 
time of recombination the universe is matter dominated r\ eq <C 77* . In this case the gravitational potentials are sourced 
by the dark matter component and their evolution is given in Sec .B 1. At linear order the gravitational potentials 
remain constant in time, while at second-order they are given by Eq. (B.4). In turn the gravitational potentials act 
as an external force on the CMB photons as in the equation (4.11) describing the CMB energy density evolution in 
the tightly coupled regime. 

For the regime of interest it proves convenient to use the solution of Eq. (4.11) found in (6.9). The source functions 
5a and S v are given by Eqs. (4.5) and (4.10), respectively. In particular Sa at early times - 5a (0) appearing in 
Eq. (6.9)- vanishes. For a matter dominated period 

5a (R = 0) = (A<f)' - f^dMp* + y «)' + ffo - TH)**™**® , (8-1) 
where we have used the linear evolution equations (3.1) and (3.8) with <j>W = and 
S V {R = 0) = ^vW'djvW + Wa^ 2 - 2d 1 ^ 2 - fWa'A^ + \{ti- Vi)d j *W tyv^ - 2<J' - -d 3 Ii {2 ^ . (8.2) 

O 4 o 4 

As at linear order we are evaluating all our expressions in the limit R = 3ph/4p 7 — ► 0, while retaining a non- vanishing 
and constant value for R in the expression for the photon-baryon fluid sound speed entering in the sines and cosines, 
Eq. (3.15). Using the linear solutions (3.20) and (3.22) for the energy density and velocity of photons, the source 
functions in Fourier space read 



S A {R = 0) = 



n • n \ 108 . ., .32 
k 2 c s cos(fciC s 7y) sin(fc 2 c s r?) + —k 2 c s sin(fc 2 c 5 7?) - — 

ZD o 



ki 

ki 



^(?7 - r/j)ki • k 2 ( ^cos(k 2 c s r]) - ^ 



k 2 c 3 sin(/ciC s r?)cos(fc 2 c iS ?7) 



(8.3) 



and 



Sy(R = 0) 



2ik 
.2 



3 V10 
9 

To 

9 

To 



k\ 



i- — c s —k 2 sa\(kic s r]) sa\{k 2 c s -q) + -k % 



ki 
10 



6 iv \ 18 
-cos(fc 2 c s ?7) - — 



cos(k 2 c 8 r]) 
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-2w* +i- 



c s (jy - »ji)ki • k 2 -^ sin(fc 2 c s 77) 
k 2 



k\ 



k 2 



ki^i- sin(fc 1 c s 77) sin(fc 2 c s ?7) 
ki 



(0)^(0). 



(8.4) 
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In Sy we have used the expression (4.19) for the second-order quadrupole moment Ti^ 1 ^ of the photons in the tight 



coupling limit, with the velocity 



,(!) 



Notice that, for the modes crossing the horizon at rj > r\ eq , we have 



expressed the gravitational potential during the matter dominated period in terms of the initial value on supcrhorizon 
scales deep in the radiation dominated epoch as \E , ( 1 - ) = 9\E , ( 1 )(0)/T0. 

As for the second-order gravitational potentials we have to compute the combination & 2 > + \&( 2 ) appearing in 
Eq. (6.9). The gravitational potential is given by Eq. (B.4), while $( 2 ) is given by 



$ (2) = ^(2) _ q(2) ; 

according to the relation (A. 12), where for a matter-dominated period 



Q< 2 > = 5V- 4 d l d 3 (d l ^ (1) d^ {1 '>) - |v- 2 (5 fc * (1) a fe * (1) ) . 



1.5) 



(8.6) 



We thus find 



$(2) + ^ ( 2) = 2*< 2 )(o) - x - (a k ^d k ^ - ^ V" 2 ^'^ 1 ^ v 2 - hv-^d^id^d^) 



+ -w- 2 (d k ^d k ^), 

3 



which in Fourier space reads 

$(2)+tf(2) = 2*£>(0) 



^G(ki,k 2 ,k) V 2 ~ ^F(k 1; k 2 ,k) 



io ) nV(o)nv(o) 



where the kernels of the convolutions are given by Eq. (7.4) and 

10 (k-ki)(k-k 2 ) 



G(ki,k 2 ,k) = ki-k 2 



k 2 



(8.7) 



(8.8) 



(8.9) 



(2) 

In Eq. (8.7) "Jm (0) is the initial condition for the gravitational potential fixed at some time rji > r\ eq . For the regime of 
interest it corresponds to the value of the gravitational potential on superhorizon scales during the matter-dominated 
epoch. 

We are now able to compute the integrals entering in the solution (6.9). The one involving the second-order 
gravitional potentials is straightforward to compute 



4 k 

s7 s 



J\ V > (V 2 ) +¥ 2 ^ S in[kc s ( V -v')] = -^(l-cos(fcc a »j))*W(0)- 



4 k 
3~ s 



5 1 
--2-F(ki,k 2 ,k)— (1 - cos(fcc s ?7)) 



7fc 3 c3 



G(ki,k 2 ,k) (-2 + (kc s T]) 2 + 2cos(kc s r])) 



*L?(o)nv(°) 



(8.10) 



(2) 

For the two remaining integrals, in the following we will show only the terms that in the final expression for Aq 

and the second-order velocity Vj give the dominant contributions for kr\ ^> 1, even though we have perfomed a fully 
computation. The integral over the source function 5a yields 



v 9 
di]' SACOs[kc s (rj - rj')} = - 



12 ki ■ k 2 4 1 ki ■ k 2 

T7IT + 7- ^72—72^7 sm(/c 2 c s 77) + (1 <- 2) 



where (1 <-» 2) stands by an exchange of indices. The terms that have been dropped in the expression (8.11) all 
vary in time as a cosine. However we have written the first term because, upon integration over time, it will give a 
non-negligible contribution to the velocity . For the last integral we find 



4 ik 
3 kc 



- I* dr)'S v sm[kc s (ri - »/)] = 

's Jo 



272k-k 2 + fc 2 36 1 (ki -k 2 )(k-k 2 ) 
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+ 



25 c s k 2 {k 2 -k 2 ) 
where the terms that have been dropped are proportional to cosines. 



rjsm(k 2 c s ri) + (1 <-» 2) 



^(OK^O), 
(8.12) 
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From the general solution (6.9) and the expression (B.4) for the second-order gravitational potential Vl^ 2 ) we thus 
obtain 



.(2) 
*00 



3c 2 

2 



*<?(0) + 



2(9a NL - 7)4>)*£(0) + 3^*^(0) 



cos(fcc s ?y) 



7© (i-ll)^ 1 ^ 2 '^ 2 ^^^^- 



3.13) 



We warn the reader that in writing Eq. (8.13) we have kept all those terms that contain the primordial non-Gaussianity 
parametrized by onl, and the terms which dominate at late times for krj 1. 



A. Initial conditions for the second-order gravitational potentials 

The initial condition fm (0) for the modes that cross the horizon after the equality epoch is fixed by the value of the 
gravitational potential on superhorizon scales during the matter dominated epoch. To compute this value we use the 
conservation on superhorizon scales of the curvature perturbation defined in Eq. (6.3). For a matter-dominated 
period the curvature perturbation on large-scales turns out ot be 



= -*£)(()) 



1 (5( 2 )/ 



3 Pm + 9 m 1 J ' 



(8.14) 



where we used the energy continuity equation 5^ p' m + 3TCS^p m — 3p m ^^ = and the (0 — 0) Einstein equation 
b^Pm/Pm — ~ 2\I'( 1 ) in the superhorizon limit. 

From the (0 — 0) Einstein equation on large scales S^p m /p m = — 2$( 2 ) + 4<f>( 1 ) 2 bringing 



C (2) = _!*£>(0) - | (*W(0) - *W(0)) + y*« 2 (0) . 



3.15) 



The conserved value of C (2) is parametrized as in Eq. (6.5), C (2) = 2a NL C (1)2 = (50a N L/9)* (1)2 , with C (1) = -5* (1) /3 
on large scales after the equality epoch. At second-order the two gravitational potentials in a matter dominated epoch 
differ according to Eq. (8.6) and using Eq. (8.15) we find 



^(0) = -g(a NL -l)*W 2 (0)+(A 



2V~ 4 did j {8^(0)8^(0)) - -V- 2 (d fc * (1) (0)d fc * (1) (0)) 



(8.16) 



we have expressed the gravitational potential during the matter dominated period "J'W in terms of the initial value 
on superhorizon scales after the equality epoch as \E'( 1 ) = 9\E'^ 1 ^(0)/10. In Fourier space Eq. (8.16) becomes 



27 /9y %,k a ,k) 
_(a NL -l) + 2l-l - 



<'(0)^(0) 



(8.17) 



where F is the kernel defined in Eq. (7.4). 
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We can use the explicit expression for (0) in Eq. (8.13), still keeping only the terms that contain the primordial 
non-Gaussianity parametrized by onl , and the terms which dominate at late times for ki] 3> 1 to find 



A (2) _ 



54 2 2/^9 

— («NL - 1) - g (9a NL - 19)cos(fcc s ?7) - - ( — ) G(ki, k 2 , k)?7 2 



(8.18) 



where we have also used in Eq. (8.13) c s ~ l/\/3 (except in the argument of the cosine for the reason explained in 
Sec. Ill A). 



B. Second-order photon velocity perturbation 

The second-order velocity of the photons can be obtained from Eq. (4.9) where, as usual, we drop off R 



,(2)i 



drj 



(8.19) 
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The second-order gravitional potential in matter-dominated universe can be obtained from Eqs. (8.5)-(8.6) and 
Eq. (B.4) as 



$ (2) = $P)(0) - I ^l 1 )^*' 1 ' - y V" J SiS J '(8 i * (1) 8 j $ (1) )j 77 2 -5V- 4 c^'(d J * (1) c^ (1) ) 



+ -v- 2 (a fe * (1) a fe * (1) ). 

3 



In Fourier space this becomes 

*<»> = *W(0) + 



^G(k 1 ,k 2 ,k) ?? 2 -^F(k 1 ,k 2 ,k) 



(8.20) 



5.21) 



where the kernels of the convolutions are given by Eqs. (7.4) and (8.9). The integral over $< 2 ) in Eq. (8.19) is then 
easily computed 

^ ) (0) ?? +f^G(k 1 ,k 2 ,k)7 ? 3 - A F ( kl>k2>k )^ 



rv 

/ H 

Jo 



<9 4 $ (2) = ~ik l 
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< } (0)< } (0) 



(8.22) 



where as usual the equivalence symbol means that we are evaluating a given expression in Fourier space. For the 
integral over the source function S v we use its expression in Fourier space, Eq. (8.4), and the dominant terms for 
kr] S> 1 are 



rv 

/ drfi. 
Jo 
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SJ, s | Mfc' ( Jo ) + n^-^'k^h-- -(- 



9 V k i • k 2 ki 



2Cs 3Vioy k 2 ,> C0S ( fc2C ^) 



(8.23) 

Notice that, in order to compute this integral, we must know the second-order vector metric petturbation u/. This is 
easily obtained for a matter-dominated universe from Eq. (B.5). Using Eqs. (B.7) and (B.2) one finds 

2 



aV 2 * (1) (0)<9 J * (1) (0) - <9 J V 2 * (1) (0)<9 4 * (1) (0) 



(8.24) 



giving rise to the third term in Eq. (8.23). 

Finally for the integral over A$ some caution is nedeed. Since in the final expression for v^ 1 the dominant 
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terms at late times turn out to be proportional rj, one has to use an expression for A 00 that keep track of all 
those contributions that, upon integration, scale like r\. Thus we must use the expression written in Eq. (8.13), plus 
Eq. (8.11) and Eq. (8.12), and some terms of Eq. (8.10) that have been previously neglected in Eq. (8.13). Then we 
find for krj 3> 1 



1 n 



(2) 
00 



ik l 
~4~ 



-4*£>(0),j + (2(9a NL - 7)*£> (0)<>(0) + 8tt£>(0)) 
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ik l 
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9 V G 1812ki-k 2 54 k 2 + k- (ki + k 2 )\ 
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(8.25) 



Using Eqs. (8.22), (8.23) and (8.25) we get 
.¥ 1 



,(2)i 



k 10c, 



■ (9qnl - 19) sin(fcc s r;) + i— — ki • k 2 
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k 2 



k 1 



(8.26) 



To obtain Eq. (8.26) we have also used the explicit expression (8.17) for ^m\0) and we have kept the terms depending 
on cinl parametrizing the primordial non-Gaussianity and the terms that dominate at late times for kr\ 3> 1. 
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C. Multipole moments 

The expression for the multipole moments (5.14) due to the anisotropics generated at recombination are easily 
found. The multipole moments for the source term Eq. (5.10), can be computed similarly to Eq. (7.41), and in 
addition we keep those term which are proportional to the gravitational potentials. We thus find 



~>*£m 



r' H)- 



2£+ 1 

47T 



4tt (A$ + 4$( 2 ) + 16<i>« 2 + 4*( 1 )A«) 6 eo 6 m0 + 8 i y«(k 1 )A< 1 ) (k 2 )Y 1 * m (k 1 ) ( 5, 1 



/47T 



+ i8^fv(k 1 )A?\k 2 )(-l)- m g™?- m (£ £ 0;m ^ 0) - 10 (^J «(k 1 )i;(k 2 )(-l)- m 
x ^ ir roi (ki)*i* ma (k2)enr 2 " m l + ( 8$(1) - 4 * (1) ) A ^ ± 16(* (1) v) m S« ± „W<J« + 



47T 



A. 



(2) 
2m 
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(8.27) 



where the minus sign must be used when m = and the plus sign when m = ±1. In Eq. (8.27) v m represents the 
scalar and vortical components of the velocity perturbation 



v( 2 >(k)= W ( 2 >(k)k + £ v%^±^, 



m=±l 



(8.28) 



where form an orhtonormal basis with k (and = v^). They can be easily obtained from Eq. (8.26). Moreover 
for a generic quantity /(x)v we have indicated the corresponding scalar and vortical components with (/v) m and 
their explicit expression is easily found by projecting the Fourier modes of /(x)v along the k = e 3 and (e 2 T * e i) 
directions 



(/v) m (k) = (±) 



^« (1) (ki)/(k2)y 1 * ro (ki) 



(8.29) 



where the minus sign must be used for m = — 1, +1 and the plus sign in correspondence of m = 0. 



IX. PERTURBATION MODES WITH fc > fc eq : IMPROVED ANALYTICAL SOLUTIONS 

In Sec. VII we have computed the perturbations of the CBM photons at last scattering for the modes that cross the 
horizon at rj < r\ eq under the approximation that the universe is radiation-dominated. However around the equality 
epoch, through recombination, the dark matter component will start to dominate. In this section we will account for 
its contribution to the gravitational potential and for the resulting perturbations of the photons from the equality 
epoch onwards. This leads to a more realistic and accurate abalytical solutions for the acoustic oscillations of the 
photon-baryon fluid for the modes of interest. 

The starting point is to consider the density perturbation in the dark matter component for subhorizon modes 
during the radiation dominated epoch. Its value at the equality epoch will fix the magnitude of the gravitational 
potential at r\ eq and hence the initial conditions for the subsequent evolution of the photons fluctuations during the 
matter dominated period. At linear order the procedure is standard (see, e.g, [39] and [33]), and we will use a similar 
one at second-order in the perturbations. 



A. Subhorizon evolution of CDM perturbations for r\ < rj eq 



From the energy and velocity continuity equations for CDM it is possible to isolate an evolution equation for the 
density perturbation Sd = 5pd/pd, where the subscript d stands for cold dark matter. In Ref. [22] we have obtained 
the Boltzmann equations up to second-order for CDM. The number density of CDM evolves according to [22] 

dn± + e *+*^?« + m y >d 2e *+** lfc t,$ n d + e*+*$ lfcV $ n d = . (9.1) 
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At linear order rid = rid + Strict and one recovers the usual energy continuity equation 

4 1) '+"2 < - 3* (1),< = 0, (9.2) 
with 6^ = S^pd/pd = 5^ rid/ fid- The CDM velocity at the same order of perturbation obeys [22] 

v^'d+Uv^ = . (9.3) 

Perturbing ticdm up to second-order we find 

Sf + v d f 3^' = -2(*« + - 2t,£><#> - 2t£><$ + 6*«'tf > + (4*J> - 2$«)^ . (9.4) 

The R.H.S. of this equation can be further manipulated by using the linear equation (9.2) to replace v^ % i yielding 

Sf + v d % - 3*( 2 )' = 4#>V> - 6 (V) 2 )' + (tf > 9 )' - 2t£><$ + 2*S>t,«* , (9.5) 

where we use = In Ref. [22] the evolution equation for the second-order CDM velocity perturbation has 

been already obtained 

v {2)i 'd + Hvf* + iJ + mJ + = 2*( 1 >'«< 1)< - 2v d 1)j d jVd 1}i - 4$«$«-* . (9.6) 

At linear order we can take the divergence of Eq. (9.3) and, using Eq. (9.2) to replace the velocity perturbation, wc 
obtain a differental equation for the CDM density contrast 

a(3*«'-4 1) ')]' = -aV 2 $«, (9.7) 

which can be rewritten as 

8 d 1] " +HS ( d y =S«, (9.8) 

where 

S« = 3*(!)" + m¥ iy + V 2 $ (1) . (9.9) 

When the radiation is dominating the gravitational potential is mainly due to the perturbations in the photons, and 
a{rf) oc r/. For subhorizon scales Eq. (9.8) can be solved following the procedure introduced in Ref. [40]. Using the 
Green method the general solution to Eq. (9.8) (in Fourier space) is given by 

Safari) =d+C 2 Hv)- I" dr/SW(rf)r/(hi{W)-Mtol)), 

Jo 

where the first two terms correspond to the solution of the homogeneous equation. At early times the density contrast 
is constant with 

4 1) (0) = ^A«(0) = -^«(0), (9.11) 
having used the adiabaticity condition, and thus we fix the integration constant as 

C 1 = -3*£ ) (0)/2, (9.12) 

and C 2 = 0. The gravitational potential during the radiation-dominated epoch is given by Eq. (B.16) and it starts 
to decay as a given mode enters the horizon. Therefore the source term behaves in a similar manner and this 
implies that the integrals over rf reach asymptotically a constant value. Once the mode has crossed the horizon we 
can thus write the solution as 

^ ) (k,i ) )=4W$( 1 '(0)ln[BWfo ) ], (9.13) 
where the constants and are defined as 

/>OC 

A (1) $ (1) (0) = / dri'S W (r]')r)' , (9-14) 
Jo 
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and 

^(i)$(")l n ( B (i)) ^ _£$(i)( ) - / dr/SWtf)r/hi(kr/). (9.15) 
2 Jo 

The upper limit of the integrals can be taken to infinity because the main contribution comes form when krj ~ 1 
and once the mode has entered the horizon the result will change by a very small quantity. Using Eq. (B.16) to 
compute the source function S^, and performing the integrals in Eq. (9.14) and (9.15) one finds that = -9.0 
and B^ ~ 0.62. More accurate values found in Rcf. [40] through a full numerical integration of the equations are 
AW = -9.6 and = 0.44. 

Before moving to the second-order case, a useful quantity to compute is the CDM velocity in a radiation dominated 
epoch. From Eq. (9.3) it is given by 



\ r*t#«Mrf) = -m^Ho ^-f^ , (9-16) 



where the last equality holds in Fourier space and we have used Eq. (B.16) (and the fact that a(r/) oc r] when radiation 
dominates). 

Combining Eq. (9.5) and (9.6) we get the analogue of Eq. (9.8) at second-order in perturbation theory 

5f" +H5f = S (2) , (9.17) 

where the source function is 

S (2) = 3 ^(2)» + m ^2 Y + V 2 $(2 ) _ 2d l (^) v ^) + V 2 ^ 1)2 + 2V 2 <I>« 2 + l - a ^V 1 ) - 6 (V 1 ) 2 )' + (<5< 1)2 )' 

- ^a+2*2 ) t,r)]'- (9-18) 
In fact we write Eq. (9.17) in a more convenient way as 

5f" 3*( 2 )" - s[ + H{8f 3*< 2 >' - Sl ) = s 2 , (9.19) 
where for simplicity we have introduced the two functions 

Sl = 4#>'*w - 6 (*« 2 )' + (4 1)2 )' - 2»?X1 + *W > ( 9 - 2 °) 

and 

S2 = V 2 $^ - 28^' v ( l )% ) + V 2 ^ 1)2 + 2V 2 $« 2 . (9.21) 

In this way we get an equation of the same form as (9.8) in the variable [5^ - 3^ - ftdri's^r}')] with source s 2 
on the R.H.S.. Its solution in Fourier space therefore is just as Eq. (9.10) 

rv rv 
4 2) -3* (2) -/ dr/8 1 {r/) = C 1 +C 2 Mv)- dr]'s 2 (v') rf [ln(fcr?') - ln(kr])] . (9.22) 
Jo Jo 

As we will see, Eq. (9.22) provides the generalization of the Meszaros effect at second-order in perturbation theory. 



B. Initial conditions 



In the next two sections we will compute explicitly the expression (9.22) for the second-order CDM density contrast 
on subhorizon scales during the radiation dominate era. First let us fix the constants C\ and C 2 by appealing to the 
initial conditions. At r\ — > the L.H.S. of Eq. (9.22) is constant, as one can check by using the results of Sec. VII B 
and the condition of adiabaticity at second-order (see, e.g., Rcf. [4, 41]) which relates the CDM density contrast at 
early times on superhorizon scales to the energy density fluctuations of photons by 

4 2) (0) = §A&>(0) - I (4 1} (0)) 2 = §A&>(0) - \ ($«(0)) 2 , (9.23) 
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where in the last step we have used Eq. (9.11). Therefore we can fix C 2 = and 



d =4 2) (0) -3* (2) (0). 



(9.24) 



Eq. (6.8) gives (0)— 4^( 2 ) (0) in terms of the primordial non-Gaussianity parametrized by onl, and the expression 



for \I/( 2 )(0) have been already computed in Eq. (8.17). Thus we find (in Fourier space) 



v (2) 



2(3a NL -l)+44 



f(ki,k 2 ,k) 

fc 2 



<>)*£(o) 



and from Eqs. (9.23) we derive the initial density contrast for CDM at second-order 



5< 2 >(0) = 



§(3^-l) + 33*^-§ 



*£ } (o)*£(°) 



Eq. (9.26) togheter with Eq. (8.17) allows to compute the constant C\ as 



(9.25) 



(9.26) 



Ci=4 2) (0) -3* (2) (0) = 



y(a NL - 1) + ^ 



a)/ 



(9.27) 



C. Computation of the integrals over the source functions 



We now compute the integrals over the functions S\ and s 2 appearing in Eq. (9.22). Let us first focus on the integral 

Notice that, using the linear equations (9.2) and (9.3) for the CDM density and velocity perturbations, the function 
si(i]') can be written in a more convenient way as 



and then 



£ drf Sl ( V ') = fa)) 2 - (0)) 2 + £ drf 



(9.28) 



(9.29) 



In Eq. (9.29) all the quantities are know being first-order perturbations: the linear gravitional potential \E'( 1 ) for a 
radiation dominated era is given in Eq. (B.16), the CDM velocity perturbation corresponds to Eq. (9.16) and the 
CDM density contrast is given by Eq. (9.13). Thus the integral in Eq. (9.29) reads (in Fourier space) 



f 

Jo 



drf 



-3AW kl • ka W-MW) ln(s (i) W) 
k 1 c s ri 



+ (9(k • ki) - 27k 2 )- lCsV ' ~ sin ( fclCs7 7') sm(k s c s r]') - k 2 c s r)'cos(k 2 c s r)') 



kfc^r/' 2 



k^c^rj' 3 



*i>)*£(o). 



(9.30) 



Let us recall that we are interested in the evolution of the CDM second-order density contrast on subhorizon scales 
during the radiation dominated epoch. Therefore once we compute the integrals we are interested in the limit of 
their expression for late times (k-q ^> 1). For the first contribution to Eq. (9.30) wc find that at late times it is well 
approximated by the expression 



Jo ' H c iV 2 H c s 



2:2 [ ~ {Mhcsri)} 2 + \n(B^k 2V ) \n(k lCsV ) 



(9 



We have computed also the remaining integral in Eq. (9.30), but it turns out to be negligible compared to Eq. (9.31). 
The reason is that the integrand oscillates with an amplitude decaying in time as and thus it leads just to a 
constant (the argument is the same we used at linear order to compute the integrals in Eq. (9.10)). Thus we can write 



i) ki • k 2 



2„ 2 



Ho 



2.2 ( ~ [ln(fcic s 77)] 2 + \n{B^k 2 rf) \n{k lCs r,) 



. (9.32) 
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We now compute the integrals over the function S2(t]) given in Eq. (9.21). Since at late times </> (1)2 ~ 1/r/ 4 and 
(vfrW v^p 1 ) j <~ l/ry 3 the main contribution to the integral will come from the two remaining terms, $( 2 ) and v^p 2 , 
whose amplitudes scale at late times as 1/rj 2 

S2 ~ V 2 $< 2 )+V 2 4 1)2 . (9.33) 
Two are the integrals that we have to compute 

rv 

dr t 's 2 (r]'WH k v') , (9-34) 

io 

and the one multiplying ln(fe^) 

rv 

dr)'s 2 {r)'W- ( 9 -35) 



f 

Jo 



Let us first consider the contributions from V 2 u^ 2 . The second integral is easily computed using the expression (9.16) 
for the linear CDM velocity. We find that at late times the dominant term is 



I 



'^ k ^ Hkn) ,- ■>■■ 



*Jc7( )*k;(°) (9-36) 



The first integral can be computed by making the following approximation: we split the integral between < ki] < 1 
and kr\ > 1 and for < krj < 1 we use the asymptotic expression 



^w-lifcSj *£>(<)) (7,7/ (9.37) 

while for krj > 1 we use the limit 



The the integral for < fc?7 < 1 just gives a constant, while the integral for kr\ > 1 brings the dominant contribution 
at late times being proportional to [ln(/c7y)] 2 so that we can write 



£ dr/V^V ln(fcr/) = A fc^l^ [i n(fcr?) ] 2 *W (0) *g 



(0) (fcr?>l). (9.39) 



As far as the contribution to the integrals (9.34) and (9.35) due to V 2 "^ 2 -* is concerned we have just to keep track 
of the initial condition provided by the primordial non-Gaussianity. We have verified that all the other terms give a 
negligible contribution. This is easy to understand: the integrand function on large scale is a constant while at late 
times it oscillates with decreasing amplitudes as rj~ 2 , and thus the integrals will tend asymptotically to a constant. 
We find that 

/"di/V^y ~-9$ (2) (0), (9.40) 
Jo 



and 



£ <fo/V 2 $< 2 V ln(fcr/) - (-9 + 9 7 - 9^^) $< 2 > (0) , (9.41) 



where 7 = 0.577... is the Euler constant, and $^ 2 ^(0) is given by Eq. (7.27). 
Therefore, from Eqs. (9.39), (9.36), and (9.40)-(9.41) we find that for krj > 1 

£ d V 's 2 (r ] 'WMkv') - Hkv)} = -^^Sr H^)M?(0)*£(0) + 94>( 2 >(0) (-9 + 9 7 - 9^) in(fo,)*( 2 >(0) 
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Let us collect the results of Eqs. (9.27), (9.32) and (9.42) into Eq. (9.22). We find that for kr] > 1 



4 2) (fc?7> 1) 



-3(a NL - l)Ai ln(flifcfj) + A? ln^fcu?) M^iM) 



^ [ln(fcic s 77)] 2 



+ 



In(Bifc2»j) ln(fc lCs r/)) + (1 <-» 2) 



_9_ 



0)^(0). 



(9.43) 



Notice that in Eq. (9.22) we have neglected \E , *- 2 \ which decays on subhorizon scales during the radiation dominated 
epoch (see Eq. (7.9), and we have used Eqs. (9.11) and (9.13). Eq. (9.43) represents the second-order Meszaros effect: 
the CDM density contrast on small scales (inside the horizon) slowly grows starting from the initial conditions that, 
at second-order, are set by the primordial non-Gaussianity parameter onl- As one could have guessed the primordial 
non-Gaussianity is just transferred linearly. The other terms scale in time as a logarithm squared. We stress that 
the computation of these terms allows one to derive the full transfer function for the matter perturbations at second 
order accounting for the dominant second-order corrections. In the next section we will use (9.43) to fix the initial 

(2) 

conditions for the evolution on subhorizon scales of the photons density fluctuations Aqq after the equality epoch. 



D. Computation of A$ for r\ > rj eq and modes crossing the horizon during the radiaton epoch 

(2) 

In this section we derive the energy density perturbations Aq of the photons during the matter dominated epoch, 
for the modes that cross the horizon before equality. In Sec. VIII we have already solved the problem assuming matter 
domination for modes crossing the horizon after equality. Thus it is sufficient to take the result (8.13) and replace 
the initial conditions 



^ + ^ } (0) 



+ f (l-^GOd.ka.ky*^. 



cos(kc s i]) + Bsm(kc s T]) 

(9.44) 



where we have restored the generic integration constants A and B, Vl/W is the linear gravitational potential (which is 
constant for the matter era) and (0) represents the initial condition for the second-order gravitational potential 
fixed at some time r\i > ?] eq . Eq. (9.43) allows to fix the proper initial conditions for the gravitational potentials 
on subhorizon scales (accounting for the fact that around the equality epoch they are mainly determined by the 

CDM density perturbations). At linear order this is achieved by solving the equation for 5^ which is obtained from 
Eq. (9.7) and the (0 — 0)-Einstein equation which reads (see Eq. (A. 9)) 

+ 3W 2 *« - V 2 *« = -^H 2 ^<#> + ^A«) . (9.45) 

On small scales one neglects the time derivatives of the gravitational potential in Eqs. (9.7) and (9.45) to obtain 

*r+*tf>' = §W a tf\ (9-46) 

where we have also dropped the contribution to the gravitational potential from the radiation component. The 
solution of this equation is matched to the value that 5^P has during the radiation dominated epoch on subhorizon 
scales, Eq. (9.13), and one finds that for r\ 3> i] eq on subhorizon scales the gravitational potential remains constant 
with 



We skip the details of the derivation of Eq. (9.47) since it is a standard computation that the reader can find, for 
example, in Refs. [33, 39]. Since around r\ eq the dark matter begins to dominate, an approximation to the result (9.47) 
can be simply achieved by requiring that during matter domination the gravitational potential remains constant to a 
value determined by the density contrast (9.13) at the equality epoch 

^k^^k, (9-48) 
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from Eq. (9.45) on small scales, leading to 



(kVe q ) 2 d lV " 1 



ln(Bifef?„) (i) 
(0.13fo, e ,) 2 * k W ' 



(9.49) 



where we used a{q) oc rj 2 during matter domination and Eq. (9.13) with Ai = —9.6 and B\ = 0.44. 

At second-order we follow a similar approximation. The general solution for the evolution of the the second-order 
gravitational potential VE^ 2 - 1 for rj > r\ eq is given by Eq. (B.4). We have to determine the initial conditions for those 
modes that cross the horizon during the radiation epoch. The (0 — 0)-Einstein equation reads 

3W* (2) ' + 3H 2 $ (2) - V 2 ¥ 2) - 6H 2 ($ (1) ) 2 - 12H^W - 3 (*«') 2 + d^d^ _ ^Mv 2 ^ 1 ' = 

-^ 2 (7^ 2) + 7 A oo)- ( 9 - 5 °) 



We fix the initial conditions with the matching at equality (neglecting the radiation component) 



(9.51) 



?(2)| 



where for small scales we neglected the time derivatives in Eq. (9.50). Using Eq. (9.43) to evaluate S d lVeq 
Eq. (9.47) to evaluate ^^\r]eq) we find in Fourier space 



and 



* (2) fe) 



IntBxkrjeq) /ki-k 2 \ ln(0.15fciry eg ) ln(0.15fc 2 ?? e9 ) , . , ,„ , ln( J Bifc 2 ?? ei j) 

•j(a NL — I) — — + — 4 r^- r^- + ln(i!ifci77 ei jj- 



(0.13fc77 eg ) 2 ' V k 2 7 (0.27fcifj e ,) 2 (0.27fc 2 r, e ,) 2 ' ^ ^ ,,eq ' (0.13fcrfc,) 2 



27 fc2 k l; k 2 [ln(fc^)] 2 + 3^ 22 gH^ e j 2 _ ln(fclCs7?eg) M^i o fc2^) + (1 ^ 2) 



k\k\ {kT] e q) 2 2 C 2 fc 2 



2 (0.13fc77e,j) 2 



(Q.13^) 2 



(9 



In Eq. (9.44) the initial condition (0) is given by Eq. (9.52) and "J/W is given by Eq. (9.47). The integration 

(2) 

constants can be fixed by comparing at 77 ~ ij eq the oscillating part of Eq. (9.44) to the solution Aqq obtained for 
modes crossing the horizon before equality and for r\ < rj eq , Eq. (7.15). Thus for r\ » rj eq and k » rj~ q we find that 



A$ = -WW(Ve q ) + Acos(kc sV ) - ^(ki.ka.k)^*^^)*^^,) 



>(i)/ 



(9.53) 



where 

A = 6* (2) (0) 



6 (k • ki)(k • k 2 ) [2k 1 k 2 cos(k 1 c s rj eq )coa(k 2 c s ri eq ) - 2k 1 k 2 cos(kc s r] eq ) + (k 2 + k% - k 2 ) sm{kic 8 r] eq ) sin(fc 2 c s 7? eg )] 
c\k x k 2 cos(kc s r) eq ) 



k-4 _l_ i-4 1 U4 _ Ok- 2 h 2 — Oh 2 h 2 — 9k- 2 h 2 



(9.54) 



and \E>( 2 )(0) is given in Eq. (8.17). 



CONCLUSIONS 



In this paper we have performed an analytical investigation of the second-order CMB anisotropies generated at 
recombination. In particular we have provided analytical solutions for the acoustic oscillations of the photon-baryon 
fluid in the tight coupling limit. One of the steps of this computation requires to generalize at second-order the 
Meszaros effect, describing the evolution of the CDM density perturbations on subhorizon scales. If on one hand 
we have kept track of the primordial non-Gaussian contribution, on the other the main effort has been put on the 
determination of all the additional second-order effects arising at recombination, which are a new potential source to 
the non-Gaussianity of the CMB anisotropies. They constitute the main core of the second-order radiation tranfer 
function necessary to establish the level of non-Gaussianity in the CMB. Our results give a simplified estimate of the 
non-linear dynamics at recombination and serve as a support for a numerical study of these effects which is under 
investigation [29] and which will provide a more accurate analysis. 
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APPENDIX A: EINSTEIN EQUATIONS 

In this Appendix we provide all the necessary expressions to deal with the gravitational part of the problem we are 
interested in, that is the second-order CMB anisotropics generated at recombination and the acoustic oscillations of 
the baryon-photon fluid. The first part of the Appendix contains the expressions for the metric and Einstein tensor 
perturbed up to second-order around a flat Friedmann-Robertson- Walker background, the energy mometum tensors 
for massless (photons) and massive particles (baryons and cold dark matter), and the relevant Einstein equations. 
The second part deals with the evolution equations and the solutions for the second-order gravitational potentials in 
the Poisson gauge. According to the regimes studied in Sections VII and VIII we have considered various epochs, in 
particular the radiation and the matter dominated eras. 

1. The metric tensor 

We adopt the Poisson gauge which eliminates one scalar degree of freedom from the goi component of the metric 
and one scalar and two vector degrees of freedom from gij . We will use a metric of the form 

ds 2 = a 2 (ry) [-e^dr] 2 + 2oj i dx l dr 1 + (e" 2 * + X ij )dx l dx j ] , (A.l) 

where a(r/) is the scale factor as a function of the conformal time 77, and uji and \ij are t ne vector and tensor 
peturbation modes respectively. Each metric perturbation can be expanded into a linear (first-order) and a second- 
order part, as for example, the gravitational potential $ = <&W + $( 2 )/2. However in the metric (A.l) the choice 
of the exponentials greatly helps in computing the relevant expressions, and thus we will always keep them where it 
is convenient. From Eq. (A.l) one recovers at linear order the well-known longitudinal gauge while at second-order, 
one finds ^ = <t> (2) - 2(</>«) 2 and *( 2 ) = V (2) + 2(^ (1) ) 2 where and cf)^, ^ (with 0« = and 

tpi 1 ) = VJ/W) are the first and second-order gravitational potentials in the longitudinal (Poisson) gauge adopted in 
Refs. [4, 42] as far as scalar perturbations are concerned. For the vector and tensor perturbations, we will neglect linear 
vector modes since they are not produced in standard mechanisms for the generation of cosmological perturbations (as 
inflation), and we also neglect tensor modes at linear order, since they give a negligible contribution to second-order 
perturbations. Therefore we take u>i and Xij to be second-order vector and tensor perturbations of the metric. 

2. The connection coefficients 

Let us give our definitions for the connection coefficients and their expressions for the metric (A.l). The number of 
spatial dimensions is n = 3. Greek indices (a, (3, /j, v, ....) run from to 3, while latin indices (a, b, i, j, k, ....to, n...) 
run from 1 to 3. The total spacetime metric g^ v has signature (— , +, +, +). The connection coefficients are defined as 

pa ap 

L ^~2 9 {~dxT + ~dxT ' ~dxT ) ■ {A - 2 > 

The Riemann tensor is defined as 

pa pa pa , pa pA pa pA / A o\ 

U fipv — 1 0v,n 1 Pfi,v + 1 Xp L f)v ~ 1 Ai^ 1 0n ■ 

The Ricci tensor is a contraction of the Riemann tensor 

Rpv = R pau ' (A-4) 

and in terms of the connection coefficient it is given by 

r - r) r a — r) r a 4- r a r 17 — r a r 17 (A ^ 

J^fiv — <~> a 1 M „ Ufj_ 1 ua T i- aa 1 ^ L ens L fj,a ■ \-"~ ) 
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The Ricci scalar is given by contracting the Ricci tensor 

i? = i?V (A.6) 



The Einstein tensor is defined as 

For the connection coefficients we find 

r° = w + 



G^u — R^v — -g^uR- (A. 7) 



r° 0i = ^j + noji. 



dx 

TL = J + Huj 1 + e 2 ^ 2 * 



dxi 



r<^ ( *-*, % + I x ;, + I(g-gf), 

3. Einstein tensor 

The components of the Einstein tensor are 

-2* 

G° = [m 2 - + 3(*') 2 - e 2$+2 ' 1 ' (W** - 2V 2 *)] , (A.9) 

P 2* 1 / „«\ u i 

G\ = 2— [a 4 *' + (W - *') - ^ W + 4W 2 - 2— — , (A.10) 

e - 2 <i. ^2 _ 2 — - 2*'$' - 3(*') 2 + + 2 *') + 2 *" 

u,- [ \ a / 

2* 

+ e 2 * (<9 fe $<9 fc $ + V 2 $ - V 2 *)] 8) + — (-d l $d 3 <f> - d l dj§ + d l d^ - d^d^ + d^d^ - d l ^dj$) 
* (d^ + V) - ^ (dH + V) + ^ («xf + Ixf \^ 3 ) ■ (ATI) 
Taking the traceless part of eq. (A. 11), we find 

*-$ = Q, (A.12) 



where Q is defined by 



with 



and 



V 2 Q=-P + 3A^, (A.13) 
P = P\, (A.14) 



P\ = d^d^ + ^(d t ^>d^-d 1 ^d J ^)+4TTG N a 2 e- 2q 'T i j7 

v 2 n = d^py. (A.i5) 



The trace of Eq. (ATI) gives therefore 
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e - 2 * t H 2 _ 2 _ _ 2$ /^, _ 3 ^/-)2 + 2W ( 3lf / _ q/j + 2$ « j + __ ( 2( 9 fc $<9 fe $ + d k ^d k ^ - 2d k <S>d k ^> + 2(P - 37V)) 
a 2 T\ . (A.16) 



3 

From Eq. (A. 10), we may deduce an equation for u 



2 



- ivV + (m 2 - 2^pj lu 1 = - (s) - ^f^j (2 + {H- *') a j $) - 8^G N a 2 e- 2 *T J ) . (A.17) 



4. Energy momentum tensors 

a. Energy momentum tensor for photons 
The energy momentum tensor for photons is defined as 

d 3p pP .p v 



rpfJ, 

7 y 



2 r d*p p»p v 

7= g )W?— ! ' (A ' 18) 

where g is the determinant of the metric (A.l) and / is the distribution function. We thus obtain 

T% = -p^l + A^ + ^fj , (A.19) 

T; ^ -\e^ Pl + \v™ + A«,«') + ^e'- V , (A.20) 

^ j = h (n< , + ^ (l + A« + (A.21) 
where p 7 is the background energy density of photons and 

are the quadrupole moments of the photons. 

b. Energy momentum tensor for massive particles 
The energy momentum tensor for massive particles of mass m, number density n and degrees of freedom 

where g m is the distribution function. We obtain 

T° m = ~Pm = ~Pm (l + 6$ + ^4 2) ) , (A.24) 

TL = -e^Pmvl = -e*+*p ro + + , (A.25) 

TL 3 = Pm (V^ + vlv m ^ = p m (V^ + v^v^^j . (A.26) 
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where p m is the background energy density of the massive particles and we have included the equilibrium temperature 



T 



APPENDIX B: SOLUTIONS OF EINSTEIN EQUATIONS IN VARIOUS ERAS 



1. Matter-dominated era 

During the phase in which the CDM is dominating the energy density of the Universe, a <~ i] 2 and we may use 
Eq. (A. 16) to obtain an equation for the gravitational potential at first-order in perturbation theory (for which 
$(i) = ^(i)) 

+3H$ (1) ' = 0, (B.l) 

which has two solutions $^ = constant and = H/a 2 . At the same order of perturbation theory, the CDM 
velocity can be read off from Eq. (A. 10) 

uW^-JLa^i). (b.2) 

At second-order, using Eqs. (A. 16) and (A. 13) and the fact that the first-order gravitational potential is constant, we 
find and equation for the gravitational potential at second-order 

( 2 )" + 3Wtf (2) ' = S m , 

Sm = -d k *Wd k *W + N = -d k $Wd k $W + (d^dj^A , (B.3) 

whose solution is 

with W(rj) = Wo/a 3 (a = 1) the Wronskian obtained from the corresponding homogeneous equation. In Eq. (B.4) 

"I'm (0) represents the initial condition (taken conventionally at r/ — > 0) deep in the matter-dominated phase. 
From Eq. (A. 17), we may compute the vector perturbation in the metric 



- ivV - 3H 2 ^dj (tfSWvM - y'W 1 )*) , (B.5) 
where we have made use of the fact that the vector part of the CDM velocity satisfies the relation — 



-or. 



The matter contrast 5^ satisfies the first-order continuity equation 



Going to Fourier space, this implies that 
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* ) -* , (0) + ^ 1) , (B.7) 
where 5^ (0) is the initial condition in the matter-dominated period. 

2. Universe filled by matter and a relativistic component 

We extend the results above for the case of CDM and a relativistic component whose energy density will be indicated 
with p v . At first-order in perturbation theory the trace of the (i — j)-component of Einstein equations, Eq. (A. 16), 
yields 

*< 1 >"+ 3H*W = HQW + VqW + \ , 

3 2 p T 

v2g(1) = V^^n^.. (b.8) 

2 pT V J 

From Eq. (A. 10) the linear CDM velocity can be expressed as 

= J n -zEL {d ^W + Hd^) - l^vW . (B.9) 

3 Pm 9 PT 

At second-order using Eq. (A. 16) and Eq. (A. 13) we find 

*( 2 >" + 3W (2) ' = S m 

S m = 6 + 2* (1) '$ (1) ' - i(2<9 fe $ (1) d fe $ (1) + d k V (1) d k ¥V - 2<9 fe $ (1) <9 fc * (1) ) 

9p m \9j p T Pm 27p m y vk 

+ n 2 ^^- + -V 2 Q^ + HQ®' + + *«)V 2 Q« + 2H 2 ^^A$ U , (B.10) 

Pt 2 3 3 pt 

where at second-order wc find 



\v 2 q^ = 3v- 2 a,^' 



2 



2 2 2 j 



+ Ih 2 ^v£K<£] _ ^2^„(i)2 _ 5 fc$(i) 5fc *(i) _ I(afc$(i)a fc $(i) - a^wa^W) , (b.h) 

2 px 1 nr7 ~" 



2 Pt m 2 



where v 2 = uW 8 ?^ and one has to employ the expression (B.9). 
For the second-order vector metric perturbation we find 



ivV + (4H 2 - 2— K = -{8\ - &V- 2 dj) f-2v&( 1 )'^$(D + 2H 2 ^v^ + 4H 2 ^ A$" + %U 2p -^v^ 
2 a J pt Pt 2 px 



+ 3H 2 ^6^>v^+m 2 (^ - *W)^v™ +3H 2 ^(^ - *(D) v aW +H 2 ^c 
Pt pt Pt Pt 



(B.12) 



3. Radiation-dominated era 



We consider a universe dominated by photons and masless neutrinos. The energy momentum tensor for massless 
neutrinos has the same form as that for the photons. During the phase in which radiation is dominating the energy 
density of the Universe, a <~ r\ and we may combine Eqs. (A. 9) and (A. 16) to obtain an equation for the gravitational 
potential ty") at first-order in perturbation theory 
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3 3 



V 2g(l) = ^^Wi^ (B13) 



where the total ansisotropy stress tensor is 



n^ = |n;. + |nt,. (B.i4) 



We may safely neglect the quadrupole and solve Eq. (B.13) setting u± = $^77. Then Eq. (B.13) becomes going 
to Fourier space 



This is the spherical Bessel function of order 1 with solutions u+ = ji(kr]/\/3), the spherical Bessel function, and 
U- = n\{krj/^/2>), the spherical Neumann function. The latter blows up as r\ gets very small and we discard it on the 
basis of the initial conditions. The final solution is therefore 

_ 3 $(i) (Q) sin ( k v/V^) ~ (kr)/V3)cos(kr)/V3) 
k (k V /V3) 3 { ' ' 

where ^(^(O) represents the initial condition deep in the radiation era. 

At the same order of perturbation theory, the radiation velocity can be read off from Eq. (A. 10) 



1 (ad^^Y 

v {1)l = '— (B 17) 

7 -2K- a 

At second order, combining Eqs. (A. 9), (A. 16), we find 



$(2)" + 4«*(2)' _ I V 2 * (2) = S* 7 , (B.18) 



3 

5 7 = 4 (V 1 )') 2 + 2$«'*«' + + *«)v 2 *« - -(d k $Wd k $w + a fc *(i)a**(i) _ a k ^d k ^) 

V / 3 3 

+ HQW + \v 2 Q^ + |($« + *(D)V 2 Q« , (B.19) 



2 2 V ; V 2 

^^(fWnfJ , (B.20) 



+ ^ V^^" 9W V 2 " 
whose solution is 



where W(rj) — (a(0)/a) 4 is the Wronskian, and ^kom is the solution of the homogeneous equation. 



horn. 

The equation of motion for the vector metric perturbations reads 
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- -VV + 4WV 



2* (1) >$ (1) - 2W? 



PT 



PT 



PT 
'.P l+Pv , 
PT 



PT 



(B.22) 
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